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ABSTRACT 

The  University  of  Texas  San  Antonio  (UTSA)  is  developing  engine  health  monitoring 
(EHM)  technology  that  compliments  the  ongoing  and  planned  research  within  AFRL. 
The  UTSA  team,  which  also  includes  Southwest  Research  Institute,  is  focusing  on 
key  technology  areas  including  engine  probabilistic-based  diagnostics  and 
prognostics,  probabilistic  life  prediction,  and  the  materials  science  aspects  of  sensors. 

The  program  consists  of  three  distinct  but  related  task  areas  that  span  EHM  from  a 
systems  engineering  level,  to  a  specific  damage-based  life  prediction  processor,  to  a 
durability  assessment  of  sensing  materials.  Task  1  is  a  systems  level  capstone  effort 
focused  on  the  information  management,  diagnostics  and  prognostics  of  EHM 
systems.  The  objectives  are  to  develop  Bayesian  learning  and  neural  networks  for 
learning  the  unknown  aspects  of  nonlinear  engine  systems  and  sensor  sensitivity 
analysis.  Task  2  is  focused  on  developing  a  probabilistic  fracture  mechanics  model 
and  ASIC  (application  specific  integrated  circuit)  implementation  for  efficient  on¬ 
board  and  real-time  assessment  of  the  damage  state  of  critical  engine  components. 
The  effort  is  to  develop  hardware  such  that  sophisticated  probabilistic  fracture 
mechanics  algorithms  can  be  placed  on-board  for  evaluation  of  detected  defects.  Task 
3  is  focused  on  the  development  of  much-needed  durability  models  for  thin  film 
sensors  that  are  either  in  common  use  or  likely  candidates  for  monitoring  changes  in 
engine  performance  or  detecting  and  monitoring  defects  in  fracture  critical  engine 
components. 

Project  Summary 

Task  1.  Engine,  Sensor  &  Data  Communication  Systems  Modeling 

T1.1  Background  &  Objectives 

Recent  interest  in  unmanned  aircraft  and  vehicles  has  increased  the  need  to 
accommodate  failures  of  its  components  such  as  actuators  and  sensors  within  the 
control  systems.  Since  faulty  sensor  readings  can  lead  to  undesirable  situations  such 
as  dispatch  delays,  degraded  engine  performance,  and  possibly  unsafe  operation, 
sensor  validation  is  an  essential  task  in  the  safe  operation  of  an  aircraft  engine.  A 
sensor  fault  is  a  deviation  from  the  sensor’s  expected  behavior.  Experience  has 
shown  that  even  carefully  designed  and  tested  sensors  may  encounter  such  faults. 
Consequently,  timely  detection  and  validation  of  sensor  faults  are  very  important  for 
safe  and  efficient  engine  operations. 

To  mitigate  the  risk  associated  with  engine  sensor  faults,  we  have  proposed  a  model- 
based  probabilistic  sensor  model  that  is  able  to  identify  sensor  faults  and  to  correct 
them  where  there  is  redundancy  in  the  sensor  information.  Our  current  objectives 
include  designing  a  model-based  sensor  validation  model  and  engine  modeling,  and 
development  of  a  detection  mechanism  that  is  to  discover  the  occurrences  of  sensor 
faults.  Since  more  advanced  wireless  technology  capable  of  operating  in  the  high 
temperature  engine  environment  in  not  available  at  this  moment,  as  the  preparation 
for  the  future  work,  we  have  studied  Bluetooth  wireless  technology  as  the  framework 
in  this  research.  With  the  given  information  from  multiple  remote  sensors  that  needs 
to  be  processed,  the  Bluetooth  technology  is  investigated  in  order  to  assess  the  real 
performance  and  coexistence  issues  of  wireless  sensors  in  the  presence  of  mutual 
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interference.  Analytical  lower  bounds  on  the  aggregated  throughput  offered  by 
Bluetooth  technology  are  evaluated  at  the  optimal  operating  conditions. 

T1.2  Approach  and  Accomplishments 

We  have  studied  and  designed  simulation  models  that  have  included  probabilistic 
information  of  sensors.  Since  faulty  sensor  readings  can  lead  to  an  unreliable 
situation,  such  as  dispatch  delays  and  possible  unsafe  operation,  fast  and  accurate 
sensor  validation  is  a  crucial  component  for  guaranteed  system  performance.  Built- 
in-test  (BIT)  is  an  approach  becoming  popular  in  many  new  complex  systems  and  can 
likely  be  extended  to  include  sensors.  Along  these  lines,  we  describe  a  Markov-jump 
state  space  model  approach.  We  also  construct  a  neural  network  engine  simulation 
model  to  “learn  or  capture”  the  unknown  aspects  that  occur  in  a  more  traditional 
nonlinear  state-space  model.  We  develop  several  Bayesian  solutions  to  sensor 
validation,  based  on  particle  filtering  algorithms:  Mixture  Kalman  filter  and 
Stochastic-M  algorithms,  and  expectation  maximization  (EP)  algorithm.  In  general, 
the  EP  algorithm  has  less  complexity  and  better  accuracy.  The  block  diagram  of  the 
designed  system  is  shown  in  Figure  1. 


Figure  1:  Block  diagram  of  the  designed  system 

As  for  the  progress  of  wireless  sensor  technology,  we  have  explored  actual  state  of  art 
of  sensor  technology  (hardware  and  technical  literature)  for  monitoring  an  aircraft 
engine,  and  focused  on  a  model  for  the  wireless  architecture  for  short-range  data 
communications.  In  addition,  we  derived  a  completely  analytical  framework,  which 
allows  the  determination  of  the  aggregate  throughput  offered  by  a  group  of  co-located 
Piconets,  as  well  as  the  determination  of  the  optimal  operating  conditions  while  taking 
propagation  aspects  into  account.  For  data  communication,  we  also  explore  the 
performance  of  several  Piconets  located  in  the  same  region  that  follows  an 
IEEE802.11b  standard.  An  approximated  threshold-based  approach  on  Bluetooth 
Piconets  is  proposed  to  determine  the  probability  of  packet  collision.  Assuming 
interference  limited  operating  conditions  caused  by  the  presence  of  the  interference 
due  to  the  packet  losses,  we  have  developed  result  of  receiving  maximum  throughput 
given  randomly  distributed  sensor  coverage  areas. 
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Task  2.  Development  of  a  Probabilistic  Fracture  Mechanics  Model  for 
Implementation  in  an  Application  Specific  Integrated  Circuit 

T2.1  Background  &  Objectives 

The  current  state-of-the-art  in  engine  health  monitoring  (EHM)  is  focused  on  detecting 
vibration  signals,  oil  debris,  foreign  object  damage,  and  other  engine  health  indicators, 
interpreting  the  signals  and  making  diagnostic  and  prognostics  decisions.  Sometimes  a 
wear  model  is  included  to  estimate  remaining  life.  EHM  technology  at  this  point  in  time 
attempts  to  sense  the  presence  of  a  crack  through  indirect  measurements  such  as  vibration 
abnormalities.  Some  success  is  reported  in  spin  pit  testing  but  application  in  an  on-board 
setting  is  dubious. 

There  has  been  a  tremendous  advance  in  the  science  and  computational  ability  to  predict 
the  behavior  of  cracks  (size,  cycles  to  failure)  in  the  past  50  years.  There  are  now  several 
sophisticated  computer  codes  for  crack  propagation  calculations  such  as  DARWIN®, 
NASGRO®,  and  AFGROW.  These  codes  have  many  man-years  of  investment  in 
algorithm  development,  materials  characterization  and  code  verification.  However,  at 
best  these  codes  can  be  used  in  a  design  mode  or  an  after  the  fact  analysis  given  a 
download  of  the  data  from  an  engine.  Currently,  trained  neural  networks  are  used  to 
implement  a  mechanical  or  material  degradation  in  order  to  get  real-time  performance. 
However,  neural  networks  must  be  trained  for  each  particular  situation  and  cannot  model 
the  capabilities  in  today’s  and  tomorrow’s  fracture  mechanics  codes. 

In  order  to  get  real-time  performance  and  use  the  sophisticated  fracture  mechanics 
knowledge,  we  will  implement  a  probabilistic  fracture  mechanics  algorithm  in  hardware 
to  facilitate  real-time  estimation  of  crack  behavior  and  structural  integrity.  The  objective 
of  this  research  is  to  develop  and  test  an  ASIC  (Application  Specific  Integrated  Circuit) 
that  implements  a  probabilistic  fracture  mechanics  algorithm  for  real-time  assessment  of 
structural  integrity  of  engine  components.  The  initial  focus  will  be  on  gas  turbine  engine 
disks,  however  this  technology  is  pervasive  and  therefore  could  be  extended  to  other 
engine  components,  as  well  as  other  structural  systems,  e.g.,  airframes,  off-shore 
structures. 

The  algorithm,  implemented  in  hardware,  will  provide  a  real-time  probability  of  fracture 
estimate  for  the  engine  component  (initial  focus  is  the  disk)  given  loading  inputs,  material 
properties,  initial  crack  size  estimate  and  other  inputs,  see  technical  details  below,  using 
low-cycle  fatigue  (LCF)  crack  propagation  fracture  mechanics. 


T2.2  Approach  and  Accomplishments 

A  Matlab  program  was  initiated  for  prototyping  the  computational  algorithms  with  the 
initial  focus  being  the  development  of  a  rainflow  algorithm  for  determining  and 
extracting  damaging  fatigue  stress  cycles  in  a  real-time  setting.  The  rainflow  and 
crack  growth  algorithms  were  combined  and  used  to  simulate  crack  growth  as  a 
function  of  engine  usage,  and  the  crack  growth  computed  was  compared  against 
results  from  well-established  fatigue  programs  DARWIN®  and  AFGROW. 
Subsequently,  a  probabilistic  crack  growth  algorithm  has  been  implemented  by 
incorporating  a  Monte  Carlo  algorithm  to  consider  uncertainties  in  localized  stresses 
and  crack  growth  rates. 
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On-board  implementation  has  focused  on  the  implementation  of  the  deterministic  low 
cycle  fatigue  algorithm  on  a  TI  DSP  Processor.  The  processing  board  and  sample 
results  are  shown  in  the  figure  below,  see  Figures  2  and  3.  The  input  signal  represents 
the  engine  component  RPM  values  during  a  representative  flight  are  shown  in  Figure 
3a.  The  values  are  converted  into  localized  stress  levels  by  the  on-board  algorithm, 
damaging  stress  pairs  are  determined,  and  crack  growth  is  computed.  The  cumulative 
crack  growth  is  shown  in  Figure  3b. 


Figure  3a.  Input  Signal 


Figure  3b.  Crack  Growth 


Task  3.  DURABILITY  OF  SENSING  MATERIALS 

T3.1  Background  &  Objectives 

Vehicle  health  management  (VHM)  requires  on-board  sensors  that  are  typically 
actuated  using  materials  that  exhibit  shape  memory,  piezoelectric,  or  magnetostrictive 
behavior.  In  many  cases,  these  sensing  materials  are  used  in  the  form  of  multi-layered 
structures  whose  functionalities  and  performance  characteristics  are  intimately  linked 
to  their  nanocrystalline  or  amorphous  states.  Furthermore,  the  durability  and 
functionalities  of  thin-film  materials  are  known  to  depend  on  layer  thickness  and 
likely  on  test  temperature.  Unfortunately,  these  material  properties  are  generally  not 
available  for  thin-film  materials.  Since  on-board  sensing  devices  must  be  more 
durable  than  the  components  that  are  being  monitored,  there  is  a  strong  need 
determine  the  operative  damage  mechanism  and  to  characterize  the  durability  of  thin- 
film  sensing  materials  in  the  pertinent  size  scale  and  under  conditions  that  are  relevant 
to  the  service  environment,  which  usually  involves  fatigue  or  thermomechanical 
loading.  The  initial  focus  of  this  effort  will  be  on  magnetostrictive  thin  film  sensor 
materials  because  of  the  anticipated  importance  to  VHM  and  relative  lack  of 
durability  information. 
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The  objectives  of  this  task  are:  (1)  to  determine  the  durability  of  magnetostrictive 
sensing  materials  subjected  to  mechanical,  electrical,  and  thermal  loading,  (2)  to 
evaluate  the  dependence  of  the  performance  and  durability  on  layer  thickness,  and 
service  temperature,  and  (3)  to  develop  a  durability  model  to  optimize  sensor 
durability  with  respect  to  layer  thickness  and  interfacial  strength.  The  first  year  of 
this  program  has  focused  on  a  room  temperature  assessment,  while  elevated 
temperature  assessments  are  planned  for  years  two  and  three. 

T3.2  Approach  and  Accomplishments 

An  integrated  experimental/nanomechanical  modeling  approach  is  being  utilized  to 
achieve  the  task  objectives.  Innovative  processing  techniques  are  being  used  to 
deposit  thin-film  sensor  materials  on  an  engineering  alloy  substrate.  Novel 
experimental  techniques  are  being  employed  to  characterize  the  microstructure,  yield 
strength,  fatigue,  and  sensing  (i.e.,  magnetostrictive)  properties  of  the  thin-film  sensor 
materials  at  ambient  temperature.  The  effects  of  layer  thickness,  temperature,  and 
interfacial  strength  on  the  mechanical  and  functional  properties  are  being  interrogated 
and  correlated  to  the  underlying  nanocrystalline  state  in  order  to  identify  the  pertinent 
failure  mechanisms.  This  fundamental  understanding  is  being  used  to  model  the 
failure  process  at  the  nanoscale  level.  The  initial  low  temperature  results  are 
applicable  to  fan  applications  (including  the  critical  JSF  lift  fan),  whereas  the  elevated 
temperature  results  proposed  in  the  later  years  of  the  program  would  be  applicable  to 
fracture  critical  compressor  and  turbine  components. 

The  results  indicated  that  the  interface  toughness  of  FeCo/Ti-6Al-4V  measured  by  an 
indentation  technique  is  independent  of  layer  thickness  at  ambient  temperature.  The 
critical  stress  for  interface  debonding  increased  with  decreasing  layer  thickness 
according  to  a  critical  energy  release  rate  criterion.  Ion-cleaning  improves  adhesion, 
enhance  the  interface  toughness  and  increase  the  debonding  stress  of  the  thin  films,  as 
illustrated  in  Figure  4.  The  resulting  interface  is  relatively  weak  thereby  preventing 
crack  penetration  into  Ti-6A1-4V  substrate  -  thus  interfacial  strength  can  be 
engineered  to  prevent  damage  to  the  substrate.  As  a  result,  the  FeCo  thin  films  did 
not  alter  the  fatigue  life  of  the  Ti-6A1-4V  substrate.  The  presence  of  microcracks  and 
localized  interface  debonds  in  the  thin  films  did  not  affect  the  functionality  of  the 
sensor  to  detect  strain  via  the  inverse  magneto-elastic  effect.  The  overall  durability  of 
the  films  was  optimized  to  be  in  the  range  needed  for  practical  application  of  the  film 
as  an  imbedded  sensor. 
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Figure  4.  Ion-cleaning  increased  the  interface  toughness  from  2  J/m2  to  4  J/m2 
and  increased  the  debonding  stress  to  provide  a  durable  FeCo  thin 
film.  A  comparison  of  the  measured  and  predicted  debonding 
stresses  for  FeCo  thin  films  on  Ti-6AI-4V  shows  that  the  debonding 
stress  in  FeCo  thin  film  increases  with  decreasing  film  thickness 
according  to  a  critical  elastic  energy  release  rate  criterion  shown  in 
the  inserted  equation. 

Accomplishments  of  this  task  include  (1)  measuring  the  elastic  and  hardness 
properties  of  FeCo  thin  films  by  nanoindentation,  (2)  characterizing  the  interface 
properties  and  identifying  the  fracture  mechanisms  in  FeCo/Ti-6Al-4V,  (3)  optimizing 
the  processing  condition  and  the  interface  properties  to  achieve  a  durable  FeCo  thin 
film  that  did  not  cause  a  fatigue  life  debit  to  the  Ti-6A1-4V  substrate,  (4)  a 
demonstration  of  the  functionality  of  the  FeCo  thin  films  as  a  strain  and  crack 
detection  sensor  in  the  presence  of  micorcracks  and  interface  debonds,  and  (5)  the 
formulation  of  a  strategy  for  optimizing  interfacial  strength  so  that  thin  films  can  be 
engineered  to  meet  the  durability  requirement  for  embedded  sensors  without 
damaging  the  substrate. 

We  have  discovered  that  interface  toughness  of  FeCo/Ti-6Al-4V  is  independent  of 
layer  thickness  and  found  out  that  delamination  strength  of  FeCo/Ti-6Al-4V  increases 
with  decreasing  layer  thickness.  Also,  we  found  out  that  weak  interface  prevents 
crack  penetration  into  Ti-6AI-4V  substrate  -  thus  interfacial  strength  can  be 
engineered  to  prevent  damage  to  the  substrate.  We  have  discovered  that  low-cycle 
fatigue  lives  of  Ti-6A1-4V  are  identical  with  and  without  FeCo  thin  films.  Our  major 
achievement  is  that  we  are  able  to  identified  strategy  of  optimizing  interfacial  strength 
so  that  thin  films  can  be  engineered  to  meet  the  durability  requirement  for  embedded 
sensors 
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Introduction 


The  UTSA  Engine  Health  Monitoring  project  covers  three  major  task  areas  which  are 
described  here  as:  (1)  System  modeling  of  engine  and  sensor  systems  using  neural 
networks  and  Bayesian  learning.  Data  communication  is  also  an  aspect  of  the 
modeling.  (2)  Development  of  a  probabilistic  fracture  mechanics  model  and  a  design 
for  a  hardware  implementation.  (3)  The  durability  and  development  of  sensor 
materials.  For  convenience  this  report  is  organized  under  these  three  areas  and  each 
task  area  includes  a  background  &  objectives  description,  a  progress  description  for 
each  quarter,  and  references  in  each  area. 

1.  Engine,  Sensor  &  Data  Communication  Systems  Modeling 


1.1  Background  &  Objectives 

Recent  interest  in  unmanned  aircraft  and  vehicles  has  increased  the  need  to 
accommodate  failures  of  its  components  such  as  actuators  and  sensors  within  the 
control  systems.  Since  faulty  sensor  readings  can  lead  to  undesirable  situations  such 
as  dispatch  delays,  degraded  engine  performance,  and  possibly  unsafe  operation, 
sensor  validation  is  an  essential  task  in  the  safe  operation  of  an  aircraft  engine.  A 
sensor  fault  is  a  deviation  from  the  sensor’s  expected  behavior.  Experience  has 
shown  that  even  carefully  designed  and  tested  sensors  may  encounter  such  faults. 
Consequently,  timely  detection  and  validation  of  sensor  faults  are  very  important  for 
safe  and  efficient  engine  operations. 

To  mitigate  the  risk  associated  with  engine  sensor  faults,  we  have  proposed  a  model- 
based  probabilistic  sensor  model  that  is  able  to  identify  sensor  faults  and  to  correct 
them  where  there  is  redundancy  in  the  sensor  information.  Our  current  objectives 
include  designing  a  model-based  sensor  validation  model  and  engine  modeling,  and 
development  of  a  detection  mechanism  that  is  to  discover  the  occurrences  of  sensor 
faults.  With  the  given  information  from  multiple  remote  sensors  that  needs  to  be 
processed,  the  Bluetooth  technology  is  investigated  in  order  to  assess  the  real 
performance  and  coexistence  issues  of  wireless  sensors  in  the  presence  of  mutual 
interference.  Analytical  lower  bounds  on  the  aggregated  throughput  offered  by 
Bluetooth  technology  are  evaluated  at  the  optimal  operating  conditions. 

1.2  First  Quarter  Accomplishment 

We  have  investigated  both  linear  and  non-linear  principal  component  analysis  (PCA) 
of  sensor  information  in  the  first  quarter.  We  also  studied  and  designed  simulation 
models  that  have  included  probabilistic  information  of  sensors.  As  for  the  progress  of 
wireless  sensor  technology,  we  have  explored  actual  state  of  art  of  sensor  technology 
(hardware  and  technical  literature)  for  monitoring  an  aircraft  engine,  and  focused  on  a 
model  for  the  wireless  architecture  for  short-range  data  communications.  We  also 
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studied  the  software  package  of  the  Bluetooth  system  in  order  to  be  coupled  with  a 
low  speed  data  acquisition/host  processor  board. 


1.3  Second  Quarter  Accomplishment 

We  have  investigated  in  this  quarter  the  design  of  a  sensor  validation  model  for 
detection  of  the  occurrence  of  a  sensor  fault,  i.e.,  fault  detection  and  identification  of 
fault  type,  i.e.,  fault  identification.  We  also  selected  possible  inference  algorithms  that 
could  be  used  for  the  validation  model.  In  addition,  we  derived  a  completely 
analytical  framework,  which  allows  the  determination  of  the  aggregate  throughput 
offered  by  a  group  of  co-located  Piconets,  as  well  as  the  determination  of  the  optimal 
operating  conditions  while  taking  propagation  aspects  into  account. 


1.4  Third  Quarter  Accomplishment 

We  have  investigated  the  problem  of  on-line  sensor  validation.  Since  faulty  sensor 
readings  can  lead  to  an  unreliable  situation,  such  as  dispatch  delays  and  possible 
unsafe  operation,  fast  and  accurate  sensor  validation  is  a  crucial  component  for 
guaranteed  system  performance.  Built-in-test  (BIT)  is  an  approach  becoming  popular 
in  many  new  complex  systems  and  can  likely  be  extended  to  include  sensors.  Along 
these  lines  we  describe  a  Markov-jump  state  space  model  approach.  We  also  construct 
a  neural  network  engine  simulation  model  to  “learn  or  capture”  the  unknown  aspects 
that  occur  in  a  more  traditional  nonlinear  state-space  model.  For  data  communication, 
we  also  explore  the  performance  of  several  Bluetooth  Piconets  located  in  the  same 
region  that  follows  an  IEEE802.11b  standard.  An  approximated  threshold-based 
approach  on  Bluetooth  Piconets  is  proposed  to  determine  the  probability  of  packet 
collision. 

1.5  Final  Accomplishment 

We  have  developed  an  alternative  Bayesian  solution  to  sensor  validation  based  on 
expectation  maximization  (EP)  algorithm.  This  method  can  on-line  track  sensor  state 
efficiently.  We  investigate  the  application  of  the  developed  EP  algorithm  to  a  realistic 
system  of  a  planetary  rover  of  Mars.  The  results  show  that  compared  with  Mixture 
Kalman  Filter  developed  in  the  third  quarter,  the  EP  has  less  complexity  and  better 
accuracy. 

1.5.1  On-line  Sensor  Validation  with  Expectation  Propagation  (EP) 

We  lay  out  in  this  section  the  steps  of  EP  algorithm  for  sensor  validation  on  the  hybrid 
DSSM  (1).  For  easy  composition,  we  first  explain  the  notations  in  the  algorithm  of 
EP;  <z(xt|yi:t-i)  =  A' (x«|xt|t_i,  P«|t_i)js  |i,e  predictive  density  of  X/; 

</(xt|yi;t)  =  A''(xt|xi|t.Pt|t)  is  the  posterior  of  x,; 

<j(Xf|yi:7\-»)  =  A'’(x,|x_4|7',P_,|t) 

is  the  smoothing  density  of  x,  without  the  current  observation  y q(y,  |  x,)  =  A 
is  the  likelihood  of  y,  given  x,;  q(s,  |  y i  /-i)  =  Discrele(q(s,  =  1  |  yi:/-i),  •••  ,  <?(•*/  =  M  \ 
yi  ,.i))  is  the  filtering  prior  of  s,;  <7(5, |  yi  r.  -/)  =  Discrete(q(s,  =  1  |  yi  t.  -r).  •••  ,  q(si  -  M\ 
yi-T'-t ))  is  the  smoothing  prior  of  s, ;  q(s,  |  yi  ,)  =  Discrete(q(s,  =  1  |  y  1 /),  . . .  ,  q(s,  =  M  \ 
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y i  /))  is  the  posterior  of  s,  .  ^  =  ^*»9(5*|yi:«-i).V(HJfx,|i_i,HftPt|t_iH7t -I- Cv)  js  ^ 
normalizing  constant;  K,  is  the  Kalman  gain  and  Jt  is  the  smoothing  gain. 

The  algorithm  of  EP  for  sensor  validation  can  be  summarized  in  the  following: 

For  current  *  1  :  T 
•  1.  Initial  estimate 

-  Obtain  the  predictive  density 


-  Moment  matching 


where 


+  GU 

P.N-I  =  FP,_,|t.IFT +c„ 
9(s.|yt,.i)  =  £»,_|  Pi«,_ 

rf  >,  _  - —£■ - 1 - 

xf|f  -  - ^ — - L_ 


-  Obtain  the  estimated  livelihood 


*c|r  —  l  +  Kf (y»  —  H#<X,|,_i) 

Pe,e  =  (I  -  KtH,#)Pl|<_| 

Kr  =  P»|t_iH71(H.lPt|,_|H7,  +Cv)_i 


A» 

Qiyt  I  Sr  =  0 


(Pe,*)_,Xe|r  -  (P,|,_| r**,,,., 

*r‘  =  (p,„)-'  -  (P.i-.ir1 

9(s.  =  »|yi:.)/9(»<  =  «|yi:._l) 


•  2.  Loop  until  the  maximum  number  of  iteration  is  reached  or  convergence  happens 


-  For  t  *  current  :  -1  current- window  size+ 1 
♦  Kalman  smoothing  w  hen  t  <  current 


J,  =  P.yFTP-‘1|( 

*»fr  =  Xj|e  4-  Jr(X»^i|7'  —  Fxf|f) 

P«fT  =  P*M  +  ifT  “  Pr+l|r)d»T 

♦  Delete  old  likelihood 


P-»|T 

*-,|T 

Q(*x  =>|yi:T,_f) 


♦  Incorporate  the  likehhood  and  moment  matching 


(Pt,T  -  A*) 

P-r,r(P  »|TXe|r  —  Mt) 

q(st  =  t|yi^) 

7(y«|i»  =  t) 


q(*t  =  i|yii) 

_  ?<»«-»lyi=r.-  ,)V<y,|H,>_t|r  ) 

*'K 

 ?C«lyi;r.-tW(yi|y«.P>f  pt't\t 

p.i. 
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where 


xf|,  _  x_t[r  4-  K((y»  —  H,fx_t|r) 

P  “\t  =  (I  —  K<  H,,  )P_r|r 
K,  =  P.,|rH:.(H„P.„rH:(  +C.,)-1 

*  Obtain  new  likelihood  for  the  next  refinement 


/*«  = 


A,  — 

Q(yt\*t  =  *)  = 


(P* -  (P_r|r)_lX.-»|T 

A*”1  —  (Per)-1  —  (P-efr)”1 


-  For  t  *  cuueat-wmdowsize-l  current 

♦  Get  the  predictive  densities,  which  are  the  same  as  (7X8X9) 

♦  Incorporate  the  true  likelihoods  into  posterior  distribution 


P*|«  —  (Prlc-l+A,)-1 

X,M  =  Pr|e(pf  *4-  Pjj/_  |X»  |*_| ) 
q(st  =  i\y  it)  =  q(yt\Si  =  =  t|yiT,-#) 


1.5.2  Sensor  Validation  of  a  Planetary  Rover  system  with  EP 

The  developed  EP  is  applied  to  a  realistic  system  of  a  planetary  rover  of  Mars.  The 
model  is  a  high  dimensional  hybrid  state  space  model  or  called  linear  switching 
model.  It  consists  of  one  continuous  state  x,,  one  discrete  state  s,  and  an  observation 
variable  y,.  The  continuous  state  x,  includes  six  system  parameters  that  capture  the 
dynamic  evolution  of  the  automaton,  and  the  discrete  state  s,  represents  6  different 
fault  states  or  operational  modes  of  the  system.  The  reason  that  we  use  this  model  is 
that  it  has  more  realistic  meanings  than  other  mathematical  models  for  engine  fault 
diagnosis.  Particle  filters,  such  as  risk-sensitive  PF,  variable  resolution  PF,  Rao- 
Blackwellized  PF,  have  been  mentioned  [10]  and  used  on  this  model  [11].  Our 
purpose  of  using  EP  is  to  show  the  simplicity  and  comparable  result  of  EP  with 
Particle  filters. 

The  model  is  described  as  follows 

p(s,  =  m  I  =n)  =  l(n,m) 

-  x,  =F(s,, 5, +u  +  pn, 
y,  =  Hx,  +v, 

where  T (n,  m)  is  transition  matrix  of  the  discrete  state  s,;  F(.)  denotes  the  state  matrix, 
which  depends  on  previous  and  current  discrete  states;  H  is  measurement  matrix, 
which  is  independent  of  the  discrete  and  continuous  states;  u  is  the  input;  n,  and  v,  are 
Gaussian  noise  vectors.  The  discrete  state  s,  has  six  modes  that  refer  to  flat  driving, 
rock  under  front  wheel,  rock  under  middle  wheel,  rock  under  rear  wheel,  rock 
between  front  and  middle  wheel  and  rock  between  middle  and  rear  wheel, 
respectively.  The  six  parameters  of  the  continuous  state  represent  the  rocker  and 
bogey  angles  and  their  derivatives,  rock  width  and  rock  height.  Rocker-bogey  is  a 
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suspension  system.  The  model  is  used  to  describe  the  response  of  the  suspension 
system  to  driving  over  rocks  and  other  obstacles  to  simulate  situations  where  the 
rover’s  scientific  instruments  could  collide  with  an  obstacle,  or  where  the  rover  could 
become  “high-centered”  on  a  rock.  The  only  two  outputs  are  the  rocker  and  bogey 
angles.  One  difficulty  of  this  model  is  how  to  describe  the  state  matrix  F(.,.).  It  is  hard 
to  be  expressed  by  a  single  matrix,  since  it  is  a  series  of  nonlinear  sequential 
reactions,  involving  the  previous  and  current  discrete  states.  We  compare  through 
computer  simulation  the  performance  of  the  EP  with  a  particle  filter,  the  mixture 
Kalman  filter  (MKF),  developed  in  the  3rd  quarter.  In  our  experiment,  we  used  500 
consecutive  observations.  We  compared  the  performance  in  terms  of  accuracy  and 
complexity.  For  EP,  we  defined  the  convergence  as  the  KL  distance  between  q(s,)s 
from  two  consecutive  EP  iterations  less  than  le-3.  Based  on  this  criterion,  our 
simulation  indicated  that,  most  of  time,  EP  converges  within  2  iterations.  For  MKF, 
50  particles  were  used.  In  Figure  1.1,  we  plotted  the  number  of  fault  detection  errors 
vs.  SNR  and  it  clearly  demonstrates  the  better  performance  of  EP.  Especially,  when 
the  noise  is  large,  EP  generates  a  lot  less  errors  than  the  MKF.  We  also  compare  the 
complexity  of  the  two  algorithms.  We  considered  the  complexity  in  terms  of  the 
computer  run-time  of  for  500  time  steps  and  the  simulation  indicated  that  EP  finished 
within  12.05s  while  MKF  took  15.82s.  It  is  evident  that  EP  has  less  complexity  than 
MKF. 


Figure  1.1.  Number  of  errors  per  500  time  steps  at  different  SNRs 


1.5.3  Bluetooth  Piconets  in  the  Presence  of  Mutual  Interference 

In  previous  three  phases  of  this  project  the  coexistence  issue  has  been  investigated  in 
the  case  of  IEEE802.11b  low-rate  link  (1  Mb/s)  interfered  by  a  set  of  Bluetooth 
piconets.  We  explored  the  possibility  to  introduce  traffic  shaping  techniques,  to 
reduce  the  effect  of  interference  for  data  traffic  types.  The  effect  of  fading  was  also 
taken  into  account  by  simplified  assumptions  where  some  performance  measures, 
such  as  the  packet  error  rate,  the  number  of  packet  retransmissions  and  the 
transmission  latency,  were  introduced. 

The  goal  of  this  task  is  to  investigate  the  analytical  coexistence  of  Bluetooth  piconets 
taking  the  interference  issues  into  account  in  a  more  complicated  setting  and  to 
provide  analytical  tools  for  performance  evaluation.  Generally  speaking,  to  obtain  the 
overall  system  performance,  a  bit  level  simulation  combined  with  a  protocol  level 
simulation  does  not  represent  a  practical  solution  due  to  prohibitive  computational 
time  required.  Hence,  there  is  the  necessity  to  introduce  new  accurate  and  simple 
methods  to  performance  evaluation  that  avoid  simulation.  Based  on  the  previous 
analytical  results  in  Q1-Q3  reports  we  assume  that  only  a  terminal  can  transmit  at  a 
time  in  a  piconet,  the  number  of  7-I  piconets  surrounding  the  reference  piconet 
coincide  with  the  number  of  potentially  interfering  transmission. 

We  are  interested  in  the  average  packet  error  probability  experienced  by  the  reference 
master  which  is  given  by 


/>ep=l-(l-/>err 


Eq.  1.5.3a 


where: 

a)  Perr  is  the  probability  of  packet  loss  probability,  derived  in  Q3  report, 

b)  Ru  is  the  maximum  radius  of  a  circular  region  of  77  unsynchronized  piconets 
distributed  uniformly  in  a  circular  region  of  radius  R , 

c)  p  is  the  probability  of  frequency  collision,  derived  in  Q2  report,  and 

d)  Interferer  activity  factor  is  denoted  by  r  as  in  Q1  report. 

e)  K  is  a  properly  protection  radio  factor. 

The  results  are  plotted  in  Figure  1.2. 
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Figure  1.2.  Comparison  between  Pep  curves  for  rj  synchronized  and 
unsynchronized  piconets 

From  the  above  figure  we  find  out  that  if  the  piconets  are  synchronized  then  the 
probability  of  error  is  decreased.  For  practically  issues  and  implementation  of  sensors 
we  suggest  using  the  synchronism  of  the  transmission  in  order  to  diminish  the  mutual 
interference.  The  relation  in  Eq.  1.5.3a  gives  analytical  relation  so  we  don’t  need  to  do 
any  simulation. 

We  are  also  interested  in  the  optimization  of  the  communication  performance  of  the 
sensors  in  the  context  of  information  theory.  This  is  quantized  by  the  throughput 
bounds  Sbound  given  below  in  two  cases  when  KRU  <R,  (D  represents  the  length  of 
the  adopted  packet  type)  DM1  or  DH1,  since  we  are  considering  one-slot 
transmission). 


Eq.  1.5.3b 


and  when  KRU  >R 


In  order  to  compare  our  results  with  the  results  predicted  by  information  theory  we 
need  to  plot  these  above  results  in  the  limiting  regime  i.e  R/  Rfj  =m,  m  =  1,2,3 . 

The  results  are  plotted  in  Figure  1.3. 
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Figure  1.3.  Throughput  in  bits/sec  of  7  piconets  in  the  limiting  regime  and 
considering  mutual  interference 


The  above  figure  is  very  important  since  it  tells  us  how  many  sensors  we  can  arrange 
or  we  can  randomly  distribute  in  a  given  area  in  order  to  get  the  maximum  of  the 
throughput.  We  assumed  interference  limited  operating  conditions  caused  by  the 
presence  of  the  interference  due  to  the  packet  losses.  Including  other  effects  as  exact 
position,  the  maximum  communication  rate  for  each  sensor,  sum  capacity, 
information  rate  capacity,  user  sensors  capacity,  the  best  receiver  for  the  master  are 
open  problems  and  are  left  for  future  research. 

It  is  of  interest  to  calculate  the  exact  the  number  of  sensor  which  maximize  the  above 
throughputs.  The  results  are  given  below,  for  the  KR„  <R 


and  when  KRU  >R 
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Eq.  1.5.3d 
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Future  Research 


Eq.  1.5.3e 


Our  models  and  analytical  results  developed  in  this  task  and  overall  in  this  research 
can  be  expanded  further  for  specifically  conditions,  which  take  all  the  transmission 
aspects  into  account  such  as 

•  Propagation  effects 

•  Interference 

•  Modulation  format 

•  Coding  techniques 
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And  medium  access  control  aspects 

•  Frequency  hopping 

•  Packet  structure 

•  Traffic  load 
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2.  Development  of  a  Probabilistic  Fracture  Mechanics  Model 
for  Implementation  in  an  Application  Specific  Integrated 
Circuit 

2.1  Background  &  Objectives 

The  current  state-of-the-art  in  engine  health  monitoring  (EHM)  is  focused  on 
detecting  vibration  signals,  oil  debris,  foreign  object  damage,  and  other  engine  health 
indicators,  interpreting  the  signals  and  making  diagnostic  and  prognostics  decisions. 
Sometimes  a  wear  model  is  included  to  estimate  remaining  life. 
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EHM  technology  at  this  point  in  time  attempts  to  sense  the  presence  of  a  crack 
through  indirect  measurements  such  as  vibration  abnormalities.  Some  success  is 
reported  in  spin  pit  testing  but  application  in  an  on-board  setting  is  dubious. 

Other  research  however  is  moving  towards  direct  interrogation  of  the  component 
material  to  sense  damage  that  could  pose  a  threat.  For  example,  a  new  DARPA 
initiative  on  Material  Prognosis  is  planning  to  address  the  detection  of  fatigue  damage 
or  cracking  on  a  turbine  disk  by  depositing  a  thin-film  magnetostrictive  material  on 
the  disk,  exciting  the  sensor,  creating  a  mechanical  elastic  wave  in  the  disk,  and 
recording  the  wave  reflections.  Processing  of  the  reflected  waves  can  indicate  the 
presence  of  a  crack.  If  a  defect  is  detected,  algorithms  are  needed  to  assess  the 
severity  of  the  defect. 

There  has  been  a  tremendous  advance  in  the  science  and  computational  ability  to 
predict  the  behavior  of  cracks  (size,  cycles  to  failure)  in  the  past  50  years.  There  are 
now  several  sophisticated  computer  codes  for  crack  propagation  calculations  such  as 
DARWIN®,  NASGRO®,  and  AFGROW.  These  codes  have  many  man-years  of 
investment  in  algorithm  development,  materials  characterization  and  code 
verification.  However,  at  best  these  codes  can  be  used  in  a  design  mode  or  an  after  the 
fact  analysis  given  a  download  of  the  data  from  an  engine.  Currently,  trained  neural 
networks  are  used  to  implement  a  mechanical  or  material  degradation  in  order  to  get 
real-time  performance.  However,  neural  networks  must  be  trained  for  each  particular 
situation  and  cannot  model  the  capabilities  in  today’s  and  tomorrow’s  fracture 
mechanics  codes. 

In  order  to  get  real-time  performance  and  use  the  sophisticated  fracture  mechanics 
knowledge,  we  will  implement  a  probabilistic  fracture  mechanics  algorithm  in 
hardware  to  facilitate  real-time  estimation  of  crack  behavior  and  structural 
integrity. 

Since  it  is  unlikely  that  the  complete  level  of  sophistication  of  the  today’s  fracture 
mechanics  codes  with  many  man-years  of  investment  can  be  programmed  in  hardware 
in  a  short  amount  of  time,  a  particular  challenge  will  be  to  determine  the  most 
significant  elements  of  the  probabilistic  fracture  mechanics  algorithm.  These  elements 
will  then  be  the  focus  of  the  hardware  implementation. 

The  objective  of  this  research  is  to  develop  and  test  an  ASIC  (Application  Specific 
Integrated  Circuit)  that  implements  a  probabilistic  fracture  mechanics  algorithm  for 
real-time  assessment  of  structural  integrity  of  engine  components.  The  initial  focus 
will  be  on  gas  turbine  engine  disks,  however  this  technology  is  pervasive  and 
therefore  could  be  extended  to  other  engine  components,  as  well  as  other  structural 
systems,  e.g.,  airframes,  off-shore  structures. 

The  algorithm,  implemented  in  hardware,  will  provide  a  real-time  probability  of 
fracture  estimate  for  the  engine  component  (initial  focus  is  the  disk)  given  loading 
inputs,  material  properties,  initial  crack  size  estimate  and  other  inputs,  see  technical 
details  below,  using  LCF  crack  propagation  fracture  mechanics. 

The  objectives  are  phased  over  three  years  with  the  effort  of  years  2  and  3  dependent 
on  the  outcome  of  previous  year’s  efforts.  A  flowchart  of  the  deterministic  fracture 
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mechanics  is  shown  in  Figure  2.1.  The  elements  of  the  algorithm  are  discussed  as 
follows: 

A)  Input  from  existing  sensors  (RPM,  temperature,  vibration,  power  level 

angle) 

B)  Crack  sensor  input  -  estimates  location  and  size 

C)  Stress  transformation  -  determines  stresses  from  sensor  inputs  from  A. 

D)  Rainflow  analysis  —  determines  damage  increments 

E)  Material  property  inputs  -  crack  growth  data  (da/dN) 

F)  Fatigue  life  model  -  calculate  crack  growth  increment 

G)  Input  from  crack  sensor  regarding  current  crack  size 

H)  Failure  analysis 

I)  Life  estimate 

Rainflow  is  the  name  given  to  the  method  of  determining  stress  reversals  that  will 
cause  fatigue  crack  growth.  (A  stress  reversal  represents  a  fatigue  cycle  from  a  min- 
max-min  or  max-min-max  stress  values.)  A  balanced  rainflow  algorithm,  neither 
including  nor  excluding  too  many  stress  reversals,  is  critical  to  having  a  real-time  on¬ 
board  system.  Including  too  many  stress  reversals  would  negatively  affect  the 
performance  of  the  system  and  may  prevent  a  real-time  application.  Excluding 
damaging  stress  reversals  will  be  nonconservative. 


/ 


Crack  Sensora.c 
a:  crack  length 
c:  crack  width 


Fatigue  Life  Model: 

Fatigue  Crack 

Growth  Analysis 

Material  Properties  from  baseline 
Crack-growth  data  (KJ 


Failure  Analysis 
K>Kcora  >  critical  size 
K:  stress  intensity  factor 
Kc:  fracture  toughness 


Life  Prediction 


Figure  2.1  Flowchart  of  deterministic  fatigue  algorithm 


In  the  following  sub-sections,  the  algorithms  for  data  filtering  and  rainflow  counting 
are  explained. 
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2.2  Data  Filtering  Process 


The  beginning  step  in  a  fatigue  analysis  for  variable-amplitude  loading  such  as  engine 
rpm  data  is  to  determine  a  compact  form  of  loading  data  that  is  composed  of  multiple 
sequences  by  peaks  and  valleys.  A  counting  of  rainflow  pairs  is  easy  by  using  filtered 
loading  data,  and  calculation  time  is  also  saved  in  a  real-time  application.  In  this  sub 
section,  a  new  simple  algorithm  for  data  filtering  process  is  suggested  and 
demonstrated  using  a  sample  data  and  a  real  engine  data. 

2.2.1  Procedure 

The  fundamental  idea  of  a  data  filtering  process  is  to  select  the  peaks  and  valleys  from 
a  variable-amplitude  loading.  In  order  to  accomplish  a  data  filtering  process, 
following  two  steps  are  required. 

A.  First,  eliminate  all  intermediate  values  between  peaks  and  valleys,  and 

B.  Secondly,  eliminate  all  repeating  magnitude. 

Some  technical  computer  programming  steps  are  explained  in  the  next  section  2.2.2, 
by  using  a  sample  test  data. 

2.2.2  Example 

A  sample  test  data  is  shown  in  Table  2.1  that  contains  15  loading  steps  with  positive 
stress  data. 


Table  2.1  Sample  data  of  variable-amplitude  loading 

Loading  steps 

1  2  3  4  5  6  7  8  9  10  11  12  13  14  15 

Stress  data 

022255596674489 

The  peaks  and  valleys  stress  data  can  be  obtained,  0-9-6-7-4-9  by  insight  from  the 
sample  data  in  Table  2.1.  To  obtain  these  peaks  and  valleys  numerically,  it  is  required 
to  calculate  the  numerical  difference  in  magnitude  between  two  consecutive  points.  A 
detailed  procedure  is  explained  as  followong,  and  summarized  the  results  in  Table 
2.2.  The  Diffl  and  Diff_2  in  Table  2.2,  represent  differences,  say  current  and 
previous  difference,  respectively.  The  next  step  is  to  determine  a  slope  of  each  Diff  l 
and  Diff_2  for  checking  a  slope  change  occurrence.  The  sign  1  means  a  positive 
slope,  0  for  a  zero  slope,  and  -1  for  a  negative  slope,  respectively.  The  last  step  is  to 
modify  the  slope  to  eliminate  same  intermediate  slopes,  using  the  sign  of  Diff  l  and 
Diff_2.  The  rules  for  a  modification  slope  in  the  last  step  are: 

[sign(diff_l)=l  and  sign(diff_2)=0]  1, 

[sign(diff_l)=-l  and  sign(diff_2)=0]  -1, 

[sign(diff_l)=0  and  sign(diff_2)=0]  and  [sign(diff_l-l)=l]  ^  1, 
[sign(diff_l)=0  and  sign(diff_2)=0]  and  [sign(diff_2-l)=l]  1, 
[sign(diff_l)=0  and  sign(diff_2)=0]  and  [sign(diff_l-l)=-l]  -1, 
[sign(diff_l)=0  and  sign(diff_2)=0]  and  [sign(diff_2-l)=-l]  -1, 
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The  last  raw,  End  sign  in  Table  2.2  shows  the  modified  slope.  The  peaks  and  valleys, 
0-9-6-7-4-9  can  be  obtained  by  a  selection  of  the  stress  when  a  slope  change  has 
occurred,  with  keeping  beginning  and  end  stress  values. 

Table  2.2  Data  filtering  scheme  of  variable-amplitude  loading  (15  points) 


Loading 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

Stress 

0 

2 

2 

2 

5 

5 

;  5 

9 

6 

6 

7 

4 

4 

8 

9 

DifT  1 

* 

2 

0 

0 

3 

0 

0 

4 

-3 

0 

1 

-3 

ol 

4 

1 

Diff  2 

2 

0 

0 

3 

0 

0 

4 

-3 

0 

1 

-3 

0 

4 

1 

* 

Sign(diff  1) 

* 

1 

0 

0 

1 

0 

0 

1 

-1 

0 

1 

-1 

0 

1 

1 

Sign(diff  2) 

1 

0 

0 

1 

0 

0 

1 

-1 

0 

1 

-1 

0 

1 

1 

* 

End  sign 

* 

1 

1 

1 

1 

1 

1 

-1 

-1 

1 

-1 

-1 

1 

1 

* 

points 


*  beginning  and  ending 


This  technique  is  very  simple  and  effective  to  determine  only  the  peaks  and  valleys 
for  any  kinds  of  data  sequence,  and  easy  to  write  a  program.  Figure  2.2  shows  the  data 
reduction  results  schematically  using  Matlab,  by  using  the  sample  test  data  of  Table 


Figure  2.2  Comparison  of  raw  data  filtered  data  using  15  test  data 

The  data  reduction  process  is  also  demonstrated  to  a  real  engine  data.  In  the  Figure 
2.3,  the  first  graph  represents  a  raw  stress  data  from  engine  rpm,  and  the  second  graph 
shows  only  peaks  and  valleys  data,  which  is  filtered  from  11,360  raw  data  to  256 
points.  It  can  be  understood  that  there  is  many  repeated  data  with  the  same  value  and 
with  unnecessary  intermediate  data  between  peaks  and  valleys.  The  filtered  points  are 
only  consisted  of  peaks  and  valleys  like  shown  a  zoomed  box  in  the  second  graph, 
Figure  2.3. 
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Load  steps 


Figure  2.3  Comparison  of  original  raw  data  filtered  data  using  real  engine  data 
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2.3  Rainflow  Cycle  Counting  (Range-Pair  Method) 


In  order  to  improve  rainflow  cycle  counting  for  a  real-time  data  input  and  analysis,  a 
range-pair  method  is  chosen.  The  advantages  of  range-pair  method  are  that  it  is  simple 
to  apply  an  algorithm  because  of  comparing  just  three  consecutive  points,  and  it  may 
pick  up  many  rainflow  pairs  before  sub-iteration,  to  find  remaining  rainflow  pairs. 

2.3.1  Cycle-Extraction  Method 

1 )  Define  k|,  k2,  and  k3  from  the  stress 
Acr2 

data  without  data  reordering. 

2)  Define  two  stress  ranges  using  three 
consecutive  points: 

Acr.  =K2  -0-*.|,Acr2  -  |°*3  -0-*2| 

3)  Extract  two  points  k i,  k2  (peak  and 
valley  points)  if  the  condition. 

Act,  <  Acr2  is  satisfied. 

4)  Define  three  new  consecutive  points, 

Act,  <  Act2  is  satisfied. 

5)  Define  three  new  consecutive  points  k|=k2,  k2=kt+i,  k3=kt+2  when  condition 
Act,  <  Acr2  is  not  satisfied. 

6)  Perform  until  last  points,  and  one  reversal  counting  for  the  last  three  points,  in 
other  words,  extract  k2,  k3  if  the  condition  Act,  >  Acr2  is  satisfied. 

7)  Perform  sub-iteration  using  remaining  points  which  are  called  residue,  to  find 
out  the  stress  pairs  which  are  generated  by  residue. 

2.3.2  Verification 

ASTM_  El 049-85  test  sample  data 

ASTM  test  data  is  chosen  to  validate  the  developed  range-pair  rainflow  counting 
method.  Figure  2.5  represents  the  ASTM  data  and  result,  which  has  9  original  data 
and  4  rainflow  pairs.  The  magnitude  and  cycle  of  rainflow  pairs  are  sketched  by  each 
black  triangle.  The  magnitude  of  each  rainflow  pair  is  3(A-B),  4(E-F),  8(C-D),  6(H-I) 
respectively.  Each  rainflow  pair  in  Figure  2.6  represents  UTS  A  result.  From  Figures 
2.5-2.6,  rainflow  pairs  are  exactly  matched  in  magnitude  and  cycles. 
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Figure  2.5  ASTM  rainflow  test  results  [4] 


Figure  2.6  Rainflow  test  results  using  UTSA 


Real  engine  data  ( Air-to-Ground ) 

This  sub-section  is  introduced  to  verify  the  developed  code  UTSA  to  a  real  engine 
data.  Unlike  the  sample  test  data,  there  is  a  huge  number  of  input  data.  Darwin  is  used 
to  compare  the  rainflow  pairs  with  UTSA.  The  magnitude  of  stress  pairs  is  not  equal 
in  each  program  since  different  rainflow  counting  methods  are  adopted  in  each  code 
(simplified  rainflow  method  for  Darwin,  range-pair  method  for  UTSA).  The  total 
number  of  rainflow  pairs  is  same,  132  for  both  UTSA  and  Darwin.  There  is  no  need 
that  the  total  number  of  rainflow  pairs  should  be  the  same  for  each  code.  The  big 
difference  between  UTSA  and  Darwin  rainflow  sequence  is  that  UTSA  uses  incoming 
order  of  data  by  real  time  acquisition,  but  Darwin  modifies  the  sequence  of  order  by 
swapping  the  maximum  stress  pairs  locate  at  the  beginning. 

Figure  2.7  represents  a  Act  (sigma_max -sigma  min)  comparison  from  the  rainflow 
pairs  output  between  UTSA  and  Darwin  by  ascending  order.  From  this  figure,  UTSA 
and  Darwin  show  almost  the  same  magnitude  of  sigma  diffenerce  at  each  rainflow 
pair.  The  detailed  comparisons  in  a  specific  region  are  depicted  in  Figures  2.8-2.10.  A 
sigma  difference  shows  the  exact  same  values  between  UTSA  and  Darwin  in  a  range 
of  rainflow  pairs,  1—12  (Figure  2.8),  and  UTSA  a  little  higher  than  Darwin  in  a  range 
of  rainflow  pairs,  61-72  (Figure  2.9)  and  Darwin  a  little  higher  than  UTSA  in  a  range 
of  rainflow  pairs,  121-132  (Figure  2.10). 
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Figure  2.7  Comparison  of  Act  by  ascending  order 
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Figure  2.8  Comparison  of  Act  by  ascending  order,  1-12 


Rainflow  pairs:  61-72 
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Figure  2.9  Comparison  of  Act  by  ascending  order,  61-72 
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Rainflow  pairs:  121-132 
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Figure  2.10  Comparison  of  Act  by  ascending  order,  121-132 


2.4  Fatigue  Analysis 

In  order  to  verify  the  crack  growth  by  fatigue  analysis  using  the  developed  code 
UTSA,  constant  amplitude  loading  (test  data)  and  variable-amplitude  loading  (real 
engine  data)  are  considered.  Paris  law  is  utilized  for  UTSA  to  simplify  the  fatigue 
analysis.  Darwin  and  AFGRO  are  used  to  compare  the  crack  growth  results. 

2.4.1  Paris  Law 

The  simple  algorithm  of  Paris  law  is  introduced,  and  evaluated  in  UTSA.  Stress 
intensity  is  defined  as  following. 


A K  =  Km»-Kmm  (2.4.1) 

In  a  surface  crack,  K  varies  along  the  crack  perimeter.  Therefore,  complete 
characterization  of  crack  growth  requires  description  of  both  crack  shape  (length, 
width)  and  size.  In  this  research,  a  noncircular  crack  shape  is  considered  in 
determining  K  ,  in  equation  (2.4.2). 

K  =  (2.4.2) 

Where,  /?  is  the  dimensionless  geometry  factor.  From  an  equation  (2.4.2),  Kmlxand 
Kminare  defined  as  followings. 


-  (2.4.3) 

a  .  and  a-  are  the  maximum  and  minimum  stress  at  each  rainflow  pairs  that  can  be 
obtained  from  rainflow  counting.  By  inserting  equations  (2.4.1 -2.4.3)  into  a  Paris  law, 

—  =  c_paris(\K)m-pans,  an  incremental  crack  growth  can  be  calculated. 
dN 
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Aa(c)  =  c  paris  *  A Km-par"  *  AN 
a(c)  =  a(c )  +  A  a(c) 


(2.4.4) 


The  surface  crack  shape  and  dimension  are  represented  in  Figure  2.1 1 . 


2b  (width) 


Figure  2.11  Surface  crack  shape  and  dimension  in  a  finite  plate 
2.4.2  Verification 

2.4.2. 1  Constant  Amplitude  Loading 

One  simple  example  of  constant  amplitude  loading  is  examined  to  verify  the 
developed  code  UTSA.  Geometry  dimension  and  material  properties  are  summarized 
in  Table  2.3.  The  magnitude  of  constant  loading  is  80  ksi,  with  an  initial  value  0  ksi. 


Table  2.3  Geometry  and  material  properties 


Crack  growth  rate 

Paris  law 

Initial  crack  size 

ii 

o 

II 

o 

k> 

Fracture  toughness 

58.0 

Constants  for  Paris  law 

c  paris=2.325e-09;  m  paris=3.0594 

Specimen  dimension 

width=height=10 

Crack  type 

surface  crack(Darwin:  SC02) 

Material  property 

Darwin  :  ti64  Paris  closure.mat 

Cycles 

10 

Crack  growth  rate 

Figure  2.12  represents  the  crack  growth,  with  respect  to  number  of  cycles.  From  these 
figures,  UTSA  and  Darwin  results  show  very  good  agreement.  A  normalized 
comparison  is  represented  in  Figure  2.13,  by  setting  UTSA=1. 
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Figure  2.12  Comparison  of  acrack  with  respect  to  cycles 


Figure  2.13  Normalized  comparison  of  crack  growth  with  factor, 

UTSA=1 


Kmax  (Maximum  of  stress  intensity,  (actor [ 

Figure  2.14  represents  Kmax  and  normalized  values  comparison.  UTSA  and  Darwin 
show  very  close  results. 
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Figure  2.14  Comparison  of  Kmax  and  normalized  values  with  cycle 
increase 

Dimensionless  geometry  factory  ft 

Figure  2.15  represents  a  dimensionless  geometry  factor,  p  solution.  UTSA  and 
Darwin  show  very  close  results  to  c-direction,  but  Darwin  shows  almost  constant 
value,  0.75,  for  both  directions. 


Figure  2.15  Comparison  of  p  with  flight  number  increase 
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2.4. 2.2  Variable-Amplitude  Loading  (real  engine  data:  Air  to 
Ground) 

For  a  variable-amplitude,  three  softwares  (UTSA,  Darwin,  AFGROW)  are  applied  to 
compare  the  crack  growth,  Kmax,  and  dimensionless  geometry  factor,  p.  An  initial 
crack  size  is  assumed  0.3,  and  10  flight  number  is  assumed.  Other  geometry  and 
material  properties  are  summarized  in  Table  2.4. 


Table  2.4  Geometry  and  material  properties 


Crack  growth  rate 

Paris  law 

Initial  crack  size 

a=c=0.3 

Fracture  toughness 

58.0 

Constants  for  Paris  law 

c  paris=2.325e-09;  m  paris=3.0594 

Specimen  dimension 

width=height=  1 0 

Crack  type 

surface  crack(Darwin:  SC02, 
AFGROW:Center  semi  elliptic  surface  flaw) 

Material  property 

Darwin  :  ti64  Paris  closure.mat 

AFGROW:  ti  6Ai-4V 

Flight  numbers 

10 

Crack  growth  rate 

Figures  2.16-12.18  represent  the  crack  growth  results  in  both  directions.  Darwin 
shows  faster  crack  growth,  and  AFGROW  shows  slower  crack  growth  than  others. 
UTSA  and  UTSAD  (using  Darwin  stress  pairs  and  running  UTSA  code)  show 
almost  the  same  results,  and  their  results  are  between  Darwin  and  AFGROW  results. 
There  is  a  little  weak  point  to  compare  directly  AFGROW  with  other  results  since  the 
material  used  is  not  the  same,  as  mentioned  in  Table  2.4.  By  comparing  with  Figures 
2.16-2.18,  c  crack  grows  faster  than  acrack.  A  normalized  comparison  is  represented 
in  Figure  2.17,  by  setting  UTSA=1.  Both  cracks  show  very  little  percentage 
difference,  say  within  1%. 


Figure  2.16  Comparison  of  a  crack  with  flight  number  increase 
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Figure  2.17  Comparison  of  c_crack  with  flight  number  increase 
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Figure  2.18  Normalized  comparison  of  crack  growth  with  factor,  UTSA=1 


Kmax 

Figures  2.19-2.20  represent  Kmax  and  normalized  value  comparison.  AFGROW 
shows  approximately  5%  lower  values  relative  to  other  cases,  and  Darwin,  UTSA, 
and  UTSA  D  show  very  close  results. 
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Figure  2.19  Comparison  of  Kmax  with  flight  number  increase 
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Figure  2.20  Normalized  comparison  of  Kmax  with  with  factor,  UTSA=1 

Dimensionless  geometry  factor.  B 


Figure  2.21  represents  a  dimensionless  geometry  factor,  /?  solution.  UTSA-c, 
UTSA  D-c,  Darwin-a,  and  Darwin-c  show  the  largest  value  in  the  top  line.  UTSA-a, 
UTSA_D-a,  and  AFGROW-a  show  the  smallest  value  in  the  bottom  line. 
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Figure  2.21  Comparison  of  p  with  flight  number  increase 


2.4.2.3  Various  Mission  Types 

In  this  subsection,  various  missions  (Instrument  &  Navigation,  Functional  Check 
Flight,  Live  Fire  Missies,  Targat  Tow,  Air  to  Ground,  Live  Fire,  Transition  and 
proficiency)  are  considered  about  crack  growth  rate,  Kmax  using  UTSA  and  Darwin. 
Figures  22-28  represent  a  crack  grow  rate,  and  Figures  29-35  represent  a  Kmax.  All 
resulults  show  very  good  agreement  between  Darwin  and  UTSA.  In  the  legend,  D 
means  Darwin  and  U  means  UTSA  code  results. 


Figure  2.22  Comparison  of  a  crack  growth  rate  with  flight  number  increase 
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Figure  2.23  Comparison  of  a  crack  growth  rate  with  flight  number  increase 


Figure  2.24  Comparison  of  a  crack  growth  rate  with  flight  number  increase 
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Figure  2.25  Comparison  of  a  crack  growth  rate  with  flight  number  increase 
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Figure  2.26  Comparison  of  a  crack  growth  rate  with  flight  number  increase 
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Figure  2.27  Comparison  of  a  crack  growth  rate  with  flight  number  increase 
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Figure  2.28  Comparison  of  a  crack  growth  rate  with  flight  number  increase 


33 


30 

29 

28 

27 

26 

25 

24 

23 

22 

21 

20 

Instrument  &  Navigation 

i  ■■  • — — * - i - •  »  '  * . ■* 

— ♦ —  Darw  in 

— ■ —  UTSA 

123456789  10 

Figure  2.29  Comparison  of  Kmax  with  flight  number  increase 
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Figure  2.30  Comparison  of  Kmax  with  flight  number  increase 
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Figure  2.31  Comparison  of  Kmax  with  flight  number  increase 
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Figure  2.32  Comparison  of  Kmax  with  flight  number  increase 


Figure  2.33  Comparison  of  Kmax  with  flight  number  increase 


Figure  2.34  Comparison  of  Kmax  with  flight  number  increase 
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Figure  2.35  Comparison  ofKmax  with  flight  number  increase 

2.5  Monte  Carlo  Simulation 

The  statistical  model  is  imbedded  in  a  Monte  Carlo  simulation  [6]  to  assess  the  impact 
of  uncertainty  and  variability.  In  a  deterministic  analysis,  the  loading  results  are  used 
as  an  input  to  a  damage  model,  such  as  a  fatigue  model.  The  life  of  the  engine 
therefore  can  then  be  evaluated  at  different  stress  levels.  The  use  of  Monte  Carlo 
simulation  in  the  evaluation  process  can  overcome  this  problem.  In  this  research,  the 
uncertainty  from  stress  input  is  considered,  and  the  stress  data  was  calculated  from  the 
given  rpm  data  using  following  equation  (2.5.1). 

Stress  =  ref  _  stress  *  rpm  /  ref  _  rpm  (2.5.1) 

The  chosen  reference  stress  value  and  reference  rpm  are  1 1 3ksi  and  14,300 
respectively. 

In  order  to  compensate  uncertainty  of  stress  data,  random  numbers  were  generated  by 
log  normal  distribution,  and  multiplied  by  stress  values. 

Stress  =  Stress  *  Sigma  _ error  (2.5.2) 

In  equation  (2.5.2),  a  term  of  Sigma  error  is  represented  by  standard  deviation  and 
mean  values  which  have  log  normal  distribution. 


Sid  _dev  =  -y/log(l  +  cov2 ) 

Mean  =  \og(medianJ\  +  cov2 )  -  0.5  log(l  +  cov2 ) 

The  median  and  coefficient  of  variance  are  1 .0  and  5%,  respectively. 

2.5.1  Histogram 

Figures  2.36-2.40  represent  histograms  of  a-crack,  c-crack,  crack  area  for  sample  100, 
500,  1000,  5000,  and  10000  after  10  flights.  An  initial  crack  size  was  assumed,  0.1  in. 
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Figure  2.36  Histogram  for  sample  100 
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Figure  2.37  Histogram  for  sample  500 
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Figure  2.38  Histogram  for  sample  1000 
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Figure  2.39  Histogram  for  sample  5000 
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Figure  2.40  Histogram  for  sample  10000 
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2.5.2  Kernel  Density  Distribution 

In  this  subsection,  kernel  density  distribution  is  introduced  to  plot  each  histogram  by  a 
density  distribution.  Density  distribution  is  needed  to  study  a  Bayesian  updating 
process  in  the  next  chapter.  Following  four  different  kernels  (2.5.4)  [7]  were 
compared  in  efficiency  by  defining  efficiency  of  K  as  (2.5.5). 


Epanechnikov 

Triangular 

Normal 

(2.5.4) 

Rectangular 


—  (1-—  r 2 )/ V5  for|/|<V5,  1  otherwise 
4  5 

1  - 1/|  for  |f|  <  1 ,  0  otherwise 

\  --0  ny1 

-v/2/r  6 

for  |/|  <  1 ,  0  otherwise 


eff(K)  =  {C(kr)/C(K)}5' 4 


jj={\t2K{t)dty'n{\K(tfdt}-' 


(2.5.5) 


From  comparison  of  efficiency  [7],  there  is  very  little  to  choose  between  the  various 
kernels  on  the  basis  of  mean  integrated  square  root.  In  this  research,  normal  kernel  is 
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chosen  without  specific  reason.  Figure  2.41  represents  a  data  plot  to  evaluate  data 
continuation  when  random  numbers  are  generated  for  various  samples  by  ascending 
order.  Horizontal  axis  means  the  total  sample  numbers  and  vertical  axis  means  a- 
crack  size  distribution.  Figure  2.42  shows  a  distribution  of  crack  area  using  normal 
kernel  with  respect  to  various  sample  numbers.  More  smooth  density  distribution  can 
be  obtained  as  increase  of  sample  numbers. 
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Figure  2.41  Crack  data  plot  by  ascending  order  for  various  sample  numbers. 

Distribution  of  a  density  estimation  using  an  normal  kernel  w.r.t.  different  sample  numbers 
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Figure  2.42  Distribution  of  a  density  estimation  of  crack  area  using  normal 
kernel  with  respect  to  various  sample  numbers 
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2.5.3  Test  for  Bayesian  Updating 

If  the  value  of  a  parameter  can  be  represented  as  a  continuous  PDF,  the  prior 
assumptions  on  the  parameter  can  be  formally  updated  by  Bayes’  theorem  [8].  If  e  is 
an  observed  experimental  results,  the  prior  distribution  f\9)  can  be  revised  in  the 
light  of  experimental  results  using  Bayes’  theorem.  Let  9  be  the  random  variable  for 
the  parameter  of  a  distribution,  posterior  probability  that  will  be  in  (#,,£7  +  A0)  as 


tP(e  \9,)f\0i)A9 


(2.5.6) 


where  P(e  |  £?, )  =  P(e  \  0t  <  9  <  9t  +  A9) .  In  the  limit,  this  yields 


pjemm) 

£j{e\e,)fX9,)de 


(2.5.7) 


The  term  P{s  1 0i )  is  the  likelihood  of  observing  the  experimental  outcome  e 

assuming  that  the  value  of  the  parameter  is#.  This  is  commonly  referred  to  as  the 
likelihood  function  of  #and  denoted  L(#),  and  denominator  (normalizing  constant,  k) 
is  independent  of  6 .  Then, 


f\0)  =  kL(6)f'{G) 

k=[[xmf(0)d0) 


In  order  to  test  a  Bayesian  updating  process,  two  test  cases  are  considered.  The  first 
case  is  to  use  a  textbook  sample,  and  the  second  case  is  to  apply  engine  data. 

2.5.3. 1  Test  a  Bayesian  using  a  Textbook  Sample  [8] 

Two  examples  are  considered  for  testing  a  Bayesian  updating.  The  first  case  is  to  use 
a  uniform  prior  distribution  (known  as  diffuse  prior)  and  likelihood  function,  as  seen 
in  equations  (2.5.9)  and  (2.5.10). 

f\0)  =  1.0  O<0<1  (2.5.9) 

L(0)  =  0  (2.5.10) 

By  applying  equation  (2.5.8),  the  posterior  distribution  of  0  is 

f\9)  =  29  (2.5.11) 

This  example  is  very  simple  case,  and  Figure  2.43  represents  the  result  using  a 
developed  Matlab  code. 
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Bayesian  updating 


Figure  2.43  Bayesian  updating  of  distribution  using  a  Matlab  code 

The  second  example  is  that  the  prior  distribution  function  and  the  probability  of 
detection  are  assumed  like  as  equation  (2.5.12)  and  (2.5.13),  respectively. 


/'(#)  =  208.3# 

20-125# 

0 

P(£  |0,)  =  O 
8# 

1.0 


0  <  #  <  0.06 
0.06  <#<0.16 
#>0.16 

#<0 

0<#<0.125 

#>0.125 


(2.5.12) 


(2.5.13) 


Figure  2.44  represents  a  Bayesian  updating  of  a  distribution  using  textbook  sample 
using  a  developed  Matlab  code. 


Figure  2.44  Bayesian  updating  of  distribution  of  flaw  depth  using  a  Matlab  code 
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Figure  2.45  shows  textbook  results.  By  comparing  Figures  2.44-2.45,  it  shows  good 
agreement  in  a  posterior  distribution. 


M  BAYES1AX  ATPROACII 


Figure  2.45  Bayesian  updating  of  distribution  of  flaw  depth  [8] 

2. 5. 3. 2  Test  a  Bayesian  using  a  Real  Engine  Data 

The  10lh  flight  a-crack  data  is  assumed  as  a  prior  distribution  function,  f'(8)  and  the 
10lh  flight  c-crack  data  is  assumed  as  likelihood  function  L(0) .  This  approximated 
function  of  a  probability  of  detection  is  used  for  current  research  because  of  not 
achieving  likelihood  data  from  sensor.  It  is  assumed  that  we  can  obtain  better 
posterior  distribution  as  long  as  we  can  obtain  a  real  data  from  sensor.  Figure  2.46 
represents  a  prior  and  likelihood  distribution  when  10000  samples  were  used. 
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Figure  2.46  Bayesian  updating  of  distribution  of  flaw  depth 
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A  linear  or  nonlinear  fitting  of  density  functions  is  required  to  predict  a  posterior 
distribution  since  Bayesian  updating  can  be  performed  by  multiply  a  prior  distribution 
with  a  likelihood  distribution,  as  shown  in  equation  (2.5.8).  In  order  to  multiply  by 
two  fitting  functions  within  overlapping  range  (yellow  region  in  Figure  2.46),  each 
segment  should  have  same  integral  length.  The  developed  Matlab  code  can 
automatically  recognize  the  overlapping  range  of  prior  and  likelihood  functions,  then 
calculate  posterior  distribution  using  equation  (2.5.8).  Figure  2.47  represents  a 
posterior  distribution  with  respect  to  various  sample  size. 
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Figure  2.47  Posterior  distribution  with  respect  to  various  sample  sizes 


2.6  On  Board  Testing  for  Real  Time  Demonstration 

Central  to  the  idea  of  using  embedded  processors  for  monitoring  engine  health  is  their 
ability  to  perform  heavy  computational  loads  in  real  time.  The  task  of  real  time  signal 
processing  presents  us  with  a  dilemma.  The  faster  one  collects  measurements  of  a 
continuous  real-life  process  such  as  vibrations  within  a  turbine  engine,  the  faster  a 
processor  must  perform  the  desired  function,  in  our  case  Rainflow  algorithm.  A  high 
rate  of  measurement  (or  sample  rate)  is  absolutely  necessary  to  represent  the  real  life 
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process  accurately.  It  falls  to  the  embedded  implementation  to  achieve  the  desired 
arithmetic  function  within  an  acceptable  sample  rate.  The  implementation’s  sample 
rate  performance  can  be  evaluated  by  using  code  profiling  methods.  Prior  to  the 
hardware  implementation,  concept  modeling  was  performed  on  the  algorithm  in 
Matlab  using  script  functions.  The  results  obtained  from  these  models  were  made 
available  in  numerical  and  graphical  form.  Verification  of  hardware  results  requires 
using  similar  numeric  and  graphic  methods. 

2.6.1  Methods 

Two  design  methods  were  explored.  The  first  method  was  not  used  for  the  final 
implementation.  It  involved  exploring  an  automated  conversion  from  Matlab  code  to 
C  code.  The  second  was  the  final  and  preferred  method  for  the  implementation.  It 
involved  manual  conversion  of  Matlab  code  to  C  code  for  executing  on  the  TI  DSP 
processor  TMS320C6713.  This  method  was  successful  in  meeting  the  project 
objectives;  implementation  of  fatigue  analysis  for  life  prediction  on  a  prototype  on¬ 
board  health  analysis  system.  Items  purchased  in  regard  to  the  hardware 
implementation  objective  are  listed  below,  Table  2.5. 


Table  2.5  Materials  list 


Item 

Quantity 

Cost 

Matlab  toolbox:  Embedded  Target  for  TI  DSP 

1 

$500.00 

Matlab  toolbox:  Link  for  Code  Composer 
Studio 

1 

$200.00 

TMS320C6713  DSP  Starter  Kit  (DSK)  rev  C 

1 

$400.00 

TMS320C671  DSP  Starter  Kit  (DSK) 

1 

$400.00 

Total 

$1,500.00 

2.6.2  Real  Time  Animation  Using  Code  Composer  Studio 

Rainflow  concept  was  modeled  using  Matlab  scripting.  The  results  from  Matlab  were 
used  to  verify  the  results  computed  from  the  TI  DSP  processor  which  shown  in  Figure 

2.48.  The  output  from  the  TMSC67xx  processor  was  viewed  numerically  in  the  stdout 
window  and  graphically  using  tool  under  Code  Composer  Studio,  as  seen  in  Figure 

2.49.  Figure  2.48-a  represents  an  incoming  rpm  data  that  is  equivalent  to  Figure  2.3, 
and  Figure  2.49-b  &  c  show  the  real  time  animation  plot  of  crack  growth  propagation 
to  a  and  c  directions.  From  1 1360  raw  input  data  points  sampled,  264  critical  stress 
points  were  calculated  within  the  Peak  Detection  stage.  This  is  a  97%  reduction  of 
the  raw  data.  Including  in-flight  and  post-flight  processes,  132  values  were  calculated 
for  each  parameter,  ‘a_crack’  and  ‘ccrack’.  The  parameter  ‘acrack’  increased  from 
the  initial  value  of  0.3000  to  a  final  value  of  0.3002.  The  parameter  ‘c  crack’ 
increased  from  the  initial  value  of  0.3000  to  a  final  value  of  0.3003. 
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Figure  2.48  TMS320C6713  DSP  Processor  (16  Meg  SDRAM) 


(a)  -  raw  data  (rpm)  acquisition 
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(c)  -  c  crack  propagation 

Figure  2.49  Real  time  animation  of  data  acquisition  and  crack  propagation 


Figure  2.50  represents  a  test  buffer  size.  Using  buffer  sizes  in  the  range  of  512  to  1024 
is  most  efficient.  Using  buffer  sizes  greater  than  1024  is  an  increasingly  less  efficient 
use  of  chip  memory.  These  estimates  are  based  on  one  type  of  flight  data.  Other 
benchmarks  should  be  run  to  test  different  flight  mission  types.  For  example,  a  flight 
mission  involving  high-speed  maneuvers  produces  greater  number  of  stress  reversals. 
This  translates  to  an  increase  in  data  points  that  Rainflow  must  analyze.  An  increase 
in  data  points  results  in  more  frequent  updates  to  the  crack  size  estimate.  Every 
additional  update  to  crack  size  will  consume  additional  computation  time.  Therefore, 
the  number  of  computations  possible  per  buffer  will  decrease  due  to  the  increase  in 
computation  time. 
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Figure  2.50  Computation  vs.  buffer  size 


47 


2.6.3  Sample  Rate  Performance  Analysis 

Buffer  sizes  of  256,  512,  1024,  and  2048  were  evaluated  to  determine  the  maximum 
sample  rate. 

bufT=  [256,512,1024,2048]  //  buffer  sizes 

elk  =  225e6  //  elk  speed  (Hz) 

Tclk  =  1  /  elk  //  elk  period  (s) 

C  =  [(mod  1  +  mod  2  + ....  +  mod  M)  *  buffSize]  +  ( subf  i)  *  5  +  ( subf  2)  *  10  +  ...subfn  *  X  » 


C  represents  cycles  needed  to  process  a  buffer  of  data.  To  calculate  the  best  sample 
rate  possible  the  cycles  needed  to  compute  a  buffer  of  data  are  multiplied  by  the  cycle 
period.  To  obtain  a  per-input  sample  rate,  this  buffer  sample  rate  is  divided  by  the 
size  of  the  respective  buffer. 

Tb  =  C  *  Tclk  //  max  buffer  sample  rate  (s) 

Tmax  =  Tb  /  buff  //  max  sample  rate  achievable  (s) 

The  subfl  and  subf2  sub-functions  are  only  executed  upon  the  calculation  of  an 
output.  Since  an  output  is  not  always  computed  for  every  input,  the  algorithm 
exhibits  non-linear  behavior.  The  constant  multipliers  of  5  and  10  are  derived  from 
observing  that,  for  a  given  buffer  of  flight  data,  subfl  and  subf2  were,  at  most, 
calculated  5  and  10  times,  respectively.  Since  no  two  flights  will  record  the  same 
input  data,  these  multipliers  may  not  be  the  same.  This  will  change  the  ideal  sample 
rate  from  flight  to  flight.  Therefore,  the  ideal  sample  rate  should  be  derived  from  data 
representing  a  worst  case  scenario.  In  this  case,  a  flight  with  high  engine  stresses.  A 
worst  case  sample  rate  (lowest  sample  period)  can  be  determined  by  benchmarking 
the  implementation  with  data  sets  from  different  flight  missions.  The  implementation 
achieves  a  maximum  sample  rate  of  2.14  us  when  using  a  buffer  of  2048  floating 
point  values  as  seen  in  Table  2.6. 

Table  2.6  Maximum  input  sample  rate 


Max  data  sample  rate  (us) 

256 

512 

1024 

2048 

2.17 

2.15 

2.14 

2.14 
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3.  DURABILITY  OF  SENSING  MA  TERIALS 


3.1  Background  &  Objectives 

Vehicle  health  management  (VHM)  requires  on-board  sensors  that  are  typically 
actuated  using  materials  that  exhibit  shape  memory,  piezoelectric,  or  magnetostrictive 
behavior.  In  many  cases,  these  sensing  materials  are  used  in  the  form  of  multi-layered 
structures  whose  functionalities  and  performance  characteristics  are  intimately  linked 
to  their  nanocrystalline  or  amorphous  states.  Furthermore,  the  durability  and 
functionalities  of  thin-film  materials  are  known  to  depend  on  layer  thickness  and 
likely  on  test  temperature.  Unfortunately,  these  material  properties  are  generally  not 
available  for  thin-film  materials.  Since  on-board  sensing  devices  must  be  more 
durable  than  the  components  that  are  being  monitored,  there  is  a  strong  need 
determine  the  operative  damage  mechanism  and  to  characterize  the  durability  of  thin- 
film  sensing  materials  in  the  pertinent  size  scale  and  under  conditions  that  are  relevant 
to  the  service  environment,  which  usually  involves  fatigue  or  thermomechanical 
loading.  The  initial  focus  of  this  effort  will  be  on  magnetostrictive  thin  film  sensor 
materials  because  of  the  anticipated  importance  to  VHM  and  relative  lack  of 
durability  information. 

The  objectives  of  this  task  are:  (1)  to  determine  the  durability  of  magnetostrictive 
sensing  materials  subjected  to  mechanical,  electrical,  and  thermal  loading,  (2)  to 
evaluate  the  dependence  of  the  performance  and  durability  on  layer  thickness,  and 
service  temperature,  and  (3)  to  develop  a  durability  model  to  optimize  sensor 
durability  with  respect  to  layer  thickness  and  compliant  layer.  The  first  year  of  this 
program  is  focused  on  a  room  temperature  assessment,  while  elevated  temperature 
assessments  are  planned  for  years  two  and  three. 
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3.2  First  Quarter  Accomplishment 


An  integrated  experimental/nanomechanical  modeling  approach  is  being  utilized. 
Innovative  processing  techniques  are  being  used  to  deposit  thin-film  layered  sensor 
materials  on  an  engineering  alloy  substrate  using  processing  parameters  developed  on 
a  parallel  DARPA-sponsored  program.  Novel  experimental  techniques  will 
subsequently  be  employed  to  characterize  the  microstructure,  yield  strength,  fatigue, 
and  sensing  (i.e.,  piezoelectric  or  magnetostrictive)  properties  of  the  thin-film  sensor 
materials.  These  include  nano-indentation  to  characterize  the  elastic  modulus,  yield 
strength,  fracture  toughness,  and  interface  toughness  of  selected  thin  films  as  well  as 
SEM,  TEM,  AFM,  and  optical  microscopy  to  identify  the  failure  mechanisms.  The 
effects  of  layer  thickness  and  temperature  on  the  mechanical  and  functional  properties 
will  be  interrogated  and  correlated  to  the  underlying  nanocrystalline  or  amorphous 
state  in  order  to  identify  the  pertinent  failure  mechanisms.  A  fundamental 
understanding  of  the  relationships  between  bond  coat,  layer  thickness,  and 
temperature  on  the  durability  and  functionalities  of  thin-film  sensor  materials  will  be 
developed  by  modeling  the  failure  process  at  the  nanoscale  level.  The  initial  low 
temperature  results  would  be  applicable  to  fan  and  compressor  applications,  whereas 
the  elevated  temperature  results  proposed  in  the  later  years  of  the  program  would  be 
applicable  to  turbine  discs. 


During  this  reporting  period,  efforts  were  initiated  to  develop  experimental  techniques 
for  characterizing  thin-film  sensor  materials.  Three  thin-film  magnetorestrictive 
sensor  specimens  were  provided  by  Dr.  Bruce  Lanning  at  SwRI.  Generated  as  part  of 
the  DARPA  Material  Prognosis  Program,  these  thin-film  sensor  specimens  included 
the  following: 

Sample  1 :  2  pm  Fe-Co  thin-film  on  a  Si  substrate 

Sample  2:  7  pm  Fe-Co  thin-film  on  a  Si  substrate 

Sample  3:  1  pm  Fe-Co  with  a  Cu  compliant  layer  on  a  Ti-6AI-4V  substrate 
All  films  were  fabricated  by  physical  vapor  deposition. 

The  integrity  of  the  as-deposited  thin-films  was  characterized  by  performing 
nanoindentation  tests  and  micro-scratch  tests.  The  former  was  utilized  to  determine 
the  elastic  modulus  and  hardness  of  the  thin-film  sensor  materials,  while  the  latter  was 
utilized  to  evaluate  the  bond  strength  and  fracture  toughness  of  the  thin-film/substrate 
interface.  These  tests  were  performed  using  an  MTS  nanoindentation  test 
system  (Nanolndenter®  XP)  at  MTS,  Oak  Ridge,  TN.  We  are  currently  in  the  process 
of  purchasing  such  a  test  system  for  use  in  this  project. 

The  nanohardness  tests  were  performed  using  a  Berkovich  indenter  with  MTS’ 
patented  Continuous  Stiffness  Measurement  (CSM)  technique.  With  this  technique,  a 
small  oscillation  was  superimposed  on  the  primary  loading  signal  and  the  resulting 
response  of  the  system  was  analyzed  by  means  of  a  frequency-specific  amplifier.  The 
excitation  frequency  was  held  constant  throughout  the  test  at  45  Hz.  The  excitation 
amplitude  was  controlled  such  that  the  resulting  displacement  amplitude  remains 
constant  at  2  nm.  An  array  of  10  indents  was  performed  on  each  sample;  successive 
indents  are  separated  by  approximately  75  microns.  With  the  CSM  technique,  each 
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indentation  experiment  allows  a  continuous  measure  of  the  contact  stiffness,  S.  Using 
the  dynamic  measure  of  S,  and  established  formulae  for  elastic  modulus  and  hardness, 
every  individual  indentation  experiment  yielded  elastic  modulus  and  hardness  as  a 
continuous  function  of  surface  penetration.  The  results  for  elastic  modulus  and 
hardness  vs.  penetration  for  the  ten  indentation  experiments  on  a  particular  sample 
were  averaged  using  discrete  displacement  windows  of  approximately  5  nm.  Using 
the  data  in  each  window,  an  average  and  standard  deviation  were  calculated.  Figure 
3.1  presents  preliminary  results  of  the  elastic  moduli  of  the  thin-film  specimens  as  a 
function  of  penetration  depth  into  the  thin  film.  The  low  elastic  modulus  of  Sample  3 
likely  reflects  the  thinness  of  the  Fe-Co  film,  combined  with  the  presence  of  the  Cu 
compliant  layer. 
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Figure  3.1.  Elastic  moduli  of  thin-film  sensor  materials  as  a  function  of  indent 
penetration  into  the  thin  film  surface. 


Microscratch  tests  were  also  performed  on  individual  thin-film  specimens  via 
the  MTS  nanoindenter.  During  the  scratch  tests,  load-displacement  curves  were 
obtained  for  both  axial  and  lateral  directions.  We  are  currently  in  the  process  of 
analyzing  these  results. 


3.3  Second  Quarter  Accomplishment 

Load-displacement  curves  of  microscratch  tests  performed  on  individual  thin-film 
specimens  via  the  MTS  nanoindenter  were  analyzed.  During  the  scratch  tests, 
load-displacement  curves  were  obtained  for  both  axial  and  lateral  directions.  A 
drastic  change  in  the  load-displacement  response  occurred  when  the  thin  film 
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debonded  from  the  substrate.  Thus,  the  microscratch  test  appears  to  be  a  useful  test 
for  determining  the  debonding  strength  for  toughness  for  thin  film  sensors.  However, 
more  work  is  required  to  identify  the  onset  of  debonding  from  the  load-displacement 
curve. 

The  composition  of  the  thin  film  material  was  chosen  to  be  FesoCoso  with  possible 
additions  of  V.  The  film  thickness  chosen  for  studies  were  1  and  5  pm  either  with  or 
without  a  compliant  layer.  The  dimensions  of  the  thin  film  sensors  are  0.25”  x  1”. 
Two  different  compliant  layers  are  planned,  but  the  materials  have  not  been  selected 
at  this  time. 

A  test  matrix  was  developed  for  assessing  the  durability  of  the  thin  film  material 
bonded  on  a  Ti-6A1-4V  substrate.  The  test  matrix  was  designed  on  the  basis  of  the 
expected  fatigue  life  of  Ti-6A1-4V  at  ambient  temperature  and  at  a  stress  ratio,  R, 
of  0.1,  as  shown  in  Table  1.  Since  a  disk  is  generally  designed  to  a  low-cycle 
fatigue  (LCF)  life  of  50  to  100  kilocycles,  a  stress  range  (Aa)  of  100  ksi  was  selected 
for  fatigue  testing  of  thin  film  sensor  material  according  to  the  test  matrix  shown  in 
Table  2.  Triplicate  tests  are  planned  for  Ti-6A1-4V  with  two  thicknesses  of  F050C050 
with  or  with  a  compliant  layer.  As  shown  in  Table  2,  a  total  of  25  test  specimens  are 
to  be  fabricated,  15  tests  are  planned  with  10  tests  (spares)  to  be  determined  at  a  later 
date. 

Table  1.  Predicted  Stress-Life  Behavior  for  Ti-6AI-4V  at  Various  Maximum 
Stress  (amax)  or  Stress  Ranges  (Aa) 


C^max  [ksi] 

Act  [ksi] 

R-ratio 

Expected  Life 

110 

100 

0.1 

-15,000 

100 

90 

0.1 

60,000-100,000 

95 

85.5 

0.1 

100,000+ 

Table  2.  Test  Matrix  for  Durability  Studies  of  Fe50Co5o  Thin  Film  Bounded 
to  Ti-6AI-4V 


Aa  [ksi] 

Without  Film 

With  Film,  but  no  CL 

1  pm  5  pm 

With  Film  +  CL 

CL1  CL2 

100 

XXX 

XXX 

XXX 

XXX 

XXX 

Specimens:  25  CL  =  Compliant  Layer 

Planned  Tests:  1 5 

Spares:  10  (for  additional  replication  or  exploration  of  additional  film 

thicknesses  or  compliant  layers) 


The  test  specimen  is  in  the  form  of  a  thin  strip  with  the  dimensions  shown  in 
Figure  3.2.  Three-point  bending  will  be  used  for  the  fatigue  tests.  A  schematic  of  the 
three-point  bend  fatigue  test  set-up  is  shown  in  Figure  3.3.  Also  shown  in  Figure  3.3 
are  five  locations  where  thin  film  sensors  (0.25”  x  1”)  will  be  placed  to  assess  the 
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durability  of  the  thin  film  sensors.  Three  sensors  will  be  placed  at  different  stress 
levels  near  the  center  of  the  tension  side  of  three-point  bend  specimen.  Two  sensors 
will  be  placed  at  the  free  ends  under  zero  stress  on  the  compression  side  of  the  test 
specimen.  These  two  sensors  will  be  used  to  detect  possible  crack  initiation  in  the 
sensors  and  the  substrate  at  various  stressed  locations  in  the  center  of  the  three-point 
bend  specimens.  The  performance  of  the  sensors  will  be  monitored  both  under  active 
and  passive  modes  as  a  function  of  fatigue  cycles. 

Machining  of  the  Ti-6A1-4V  test  specimens  are  currently  underway.  An  existing  test 
fixture  is  available  for  the  three-point  bend  fatigue  tests,  but  a  few  rollers  of  the 
desired  dimensions  need  to  be  fabricated.  Fabrication  of  the  rollers  is  also  underway. 
Both  the  specimens  and  the  test  fixture  are  expected  to  be  ready  for  deposition  of  thin 
film  to  commence  in  January  2005,  followed  by  fatigue  testing. 

Theoretical  modeling  of  the  durability  of  thin  film  sensor  materials  was  initiated.  The 
modeling  efforts  were  focused  on  addressing  two  issues:  (1)  the  changes  in  the 
fatigue  resistance  of  the  substrate  due  to  the  presence  of  surface-bonded  thin  film 
sensors,  and  (2)  the  resistance  to  fatigue  crack  initiation  and  growth  in  the  thin-film 
sensor  materials.  To  address  the  first  issue,  modeling  effort  was  initiated  in  the 
reporting  period  to  examine  how  a  thin  film  sensor  can  alter  the  fatigue  crack 
initiation  and  growth  lives  of  the  substrate.  The  modeling  effort  will  continue  in  the 
next  reporting  period  and  will  involve  the  use  of  a  physics-based  fatigue  crack 
initiation  and  growth  model  developed  under  AFOSR  MEANS  program  [1]  to  predict 
fatigue  crack  initiation  and  growth  lives  of  the  substrate  without  and  without  the  thin 
film  sensor.  In  the  latter  case,  the  thickness  of  the  thin  film  materials  and  the  nature  of 
residual  stresses  in  the  material  will  also  be  considered.  Theoretical  modeling  to 
address  the  second  issue  will  commence  after  the  durability  tests  of  the  thin  film 
material  are  underway. 
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Figure  3.2.  Proposed  specimen  geometry  and  dimensions  (in  inches)  for 
three-point  bend  fatigue  testing  of  FesoCo5o  thin  film  bonded  to  a 
Ti-6A1-4V  substrate.  Specimen  is  shown  without  the  thin  film 
materials. 
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Figure  3.3.  Specimen  loading  geometry  and  approximate  sensor  dimensions  and 
locations.  The  thin  film  sensors  are  0.25”  x  1”  and  are 
schematically  shown  as  rectangular  boxes  of  either  solid  lines  or 
dashed  lines.  All  dimensions  are  in  inches. 


3.4  Third  Quarter  Accomplishment 

During  this  quarter,  Fe-Co  thin-film  microstructures  and  properties  were  characterized,  and 
initial  fatigue  tests  with  these  Films  deposited  on  a  Ti-6A1-4V  alloy  typical  of  fan  and 
compressor  components  in  turbine  engines  were  performed.  Predictions  from 
micromechanical  fatigue  models  were  compared  with  the  measured  data.  This  comparison, 
combined  with  an  initial  assessment  of  the  failure  mode  of  the  films,  suggests  that  the  film 
durability  can  be  enhanced  by  improving  the  decohesion  strength  of  the  film-to-substrate 
interface.  Methods  to  improve  the  interfacial  decohesion  have  been  identified  and  will  be 
pursued  during  the  next  quarter. 

3.5  Final  Accomplishment 

Significant  progress  was  made  to  assess  the  durability  of  FeCo  thin  films  bonded  to 
Ti-6A1-4V  in  the  last  reporting  period,  which  included:  )  measuring  the  interface 
toughness,  (2)  determining  the  debonding  stresses,  )  identifying  fracture  mechanisms, 
and  (4)  characterizing  the  fatigue  life  response  of  FeCo/Ti-6Al-4V.  Detailed 
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descriptions  of  these  findings  are  presented  in  a  recent  paper  [1].  Key  findings  in  this 
investigation  are  highlighted  in  the  next  four  subsections. 

3.5.1  Interface  Fracture  in  FeCo/Ti-6AI-4V 

Indentation  interface  fracture  toughness  tests  were  performed  at  ambient  temperature 
using  a  brale  ‘C”  indenter  (120°  diamond  cone  with  a  slightly  rounded  point)  in  a 
Rockwell  hardness  tester.  The  brale  indenter  was  pushed  into  the  FeCo  thin  film  and 
the  Ti-6A1-4V  substrate  under  a  constant  applied  load.  After  indentation,  the  diameter 
of  the  indent,  the  size  of  the  debonded  region  and  the  cracked  region  on  the  thin  film 
were  measured  by  optical  microscopy  (OM)  and  scanning  electron  microscopy 
(SEM).  The  sizes  of  the  indented  and  the  cracked  region,  as  well  as  the  residual 
stresses  in  the  thin  film,  which  were  determined  separately  by  curvature 
measurements  of  the  thin  film  deposited  on  Kapton,  were  then  utilized  as  input  to  the 
interface  fracture  model  developed  by  Drory  and  Hutchinson  [2]  to  compute  the 
elastic  strain  energy  release  rate  (Gj)  at  interface  fracture  for  FeCo  thin  films  bonded 
to  Ti-6A1-4V. 

Figures  3.4(a)  and  (b)  show  the  optical  and  secondary  electron  images  of  brale  indents 
on  FeCo/Ti-6-4  and  the  corresponding  debonded  and  cracked  regions  in  the  thin  film, 
respectively.  For  both  cases,  the  largest  diameter  was  taken  to  be  the  size  of  the 
debonded  or  the  cracked  region.  A  direct  comparison  of  the  interface  toughness  using 
the  debonded  region  determined  by  OM  and  the  cracked  region  determined  by  SEM 
indicated  that  the  diameter  of  the  cracked  region  provided  interface  toughness  values 
that  were  more  reasonable  than  those  obtained  using  the  diameter  of  the  debonded 
region  measured  by  OM.  Thus,  all  interface  toughness  values  were  determined  on  the 
basis  of  the  maximum  diameter  of  the  cracked  region  measured  by  SEM,  assuming 
the  debonded  region  extended  as  far  as  the  cracked  region.  Figure  3.4(b)  shows  that 
cracks  in  the  FeCo  thin  films  are  debonded  from  the  substrate. 

Figure  3.5(a)  presents  the  results  of  the  interface  toughness  of  a  single  4  pm  thick 
FeCo  thin  film  on  Ti-6A1-4V  determined  at  various  applied  loads.  In  a  second  set  of 
experiments,  the  applied  load  was  fixed  at  15  kg,  but  the  location  was  varied  for 
indentation  measurements  along  the  length  (25.4  mm)  of  the  4  pm  thin  film,  Figure 
3.5(b).  As  shown  in  Figure  3.5(a),  most  of  the  interface  toughness  values  (Gj) 
obtained  from  the  first  set  of  experiments  were  in  the  range  of  2  -  4  J/m2,  but 
occasionally  showed  higher  values  (6-8  J/m2)  at  higher  indentation  loads  (60  -  100 
kg).  In  comparison,  the  interface  toughness  obtained  by  varying  the  indentation 
location  at  a  constant  load  (15  kg)  showed  an  interface  toughness  value  of  7.03  ±  2.41 
J/m2  compared  to  4.45  ±  1.  74  J/m2  for  those  of  30  -  100  kg  loads.  Figure  2(b).  Thus, 
the  interface  toughness  due  to  spatial  variation  was  larger  than  that  exhibited  by 
indentation  load  variation.  Since  the  interface  toughness  value  was  similar  for  most 
indentation  loads,  it  was  thought  that  the  variations  at  some  indentation  loads  were  the 
result  of  toughness  variation  with  the  indentation  location  and  not  the  effect  of  the 
indentation  load. 
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(a)  (b) 


Figure  3.4.  Optical  and  SEM  images  of  a  brale  indent  on  FeCo/Ti-6AI-4V 
showing  debonded  and  cracked  ones  in  the  indented  thin  film:  (a) 
optical  image  taken  with  side  lighting,  and  (b)  SEM  image. 
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Figure  3.5.  Interface  toughness  of  a  4pm  FeCo  thin  film  measured  using  a  brale 
indentor  in  a  Rockwell  hardness  tester:  (a)  as  a  function  of  indentation 
load,  and  (b)  as  a  function  of  location. 


In  addition  to  4  pm  thick  films,  the  interface  toughness  values  of  three  FeCo  thin 
films  along  the  length  of  the  films  were  also  measured  for  film  thickness  of  4.56, 
3.27,  and  0.9  pm.  For  these  three  cases,  the  interface  toughness  varies  with  location 
and  ranges  from  the  lowest  at  1 .9  J/m2  to  the  highest  at  5.9  J/m2.  The  locations  of  the 
lowest  and  highest  values  of  the  interface  toughness  are  different  in  the  thin  films 
examined.  On  the  other  hand,  the  mean  and  the  standard  deviations  do  not  show  large 
variations  with  film  thickness,  as  shown  in  Figure  3  which  presents  the  interface 
toughness  as  a  function  of  film  thickness.  The  results  indicate  that  the  average  value 
of  the  interface  toughness  is  approximately  the  same  (w  4  J/m2)  with  a  standard 
deviation  in  the  range  of  ±  1  —  1.8  J/m2,  Table  1.  The  corresponding  coefficient  of 
variation  (COV)  ranges  from  23.8%  to  41.8%.  Neither  the  mean  nor  the  coefficient 
of  variation  shows  an  obvious  dependence  on  the  film  thickness  for  the  range  of 
thickness  (1-5  pm)  examined.  Thus,  it  is  concluded  that  for  the  FeCo  thin  films 
examined,  the  interface  toughness  is  independent  of  the  film  thickness. 
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Figure  3.6.  Interface  toughness  of  FeCo  thin  films  as  a  function  of  film  thickness. 


Table  1.  Interface  Toughness  (G.)  and  Debonding  Stress  (ctd)  of  FeCo  Thin  Films 
on  Ti-6A1-4V  Compared  Against  Fatigue  Lives  of  Ti-6A1-4V  with  and 
without  FeCo  Thin  Films  Tested  at  a  Maximum  Stress  of  758  MPa  and  at 
R  =  0. 1 .  Standard  Deviations  Based  on  Three  Tests  or  More  are  Shown 
in  Parentheses. 


Film  Thickness,  pm 

Gj,  J/m2 

CTd,  MPa 

Fatigue  Life,  Cycles 

0.9 

3.25(1.12) 

>758 

74,646  (7156) 

3.27 

4.53  (1.08) 

>758 

— 

3.8 

— 

668.8 

— 

4.0 

4.36(1.82) 

533.4  (65.2) 
500.7  (44.5) 

60,239  (4450) 

4.2 

— 

587.2  (46.9) 

— 

4.56 

3.56(1.0) 

323.0  (58.9) 

— 

6.8 

— 

432.5  (86.8) 

— 

7.5 

— 

305.7(133.3) 

57,380  (2460) 

Without  thin  film 

— 

— 

62,894  (3680) 

Nanoindentation  tests  were  performed  at  UTSA  to  characterize  the  elastic  modulus, 
hardness,  and  yield  strength  of  FeCo  thin  films  on  Ti-6A1-4V  for  two  different  film 
thickness  and  surface  conditions.  The  nanoindentation  tests  were  carried  out  using  the 
continuous  stiffness  measurement  technique  of  MTS  with  a  Berkovich  indenter.  An 
array  of  five  indents  was  performed  on  each  specimen  and  the  results  are  summarized 
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in  Table  2.  In  addition,  a  cross-sectional  nanoindentation  technique,  which  involved 
indenting  on  the  edge  of  a  coated  specimen  to  cause  debonding  of  the  thin  film  from 
the  substrate,  was  utilized  to  measure  the  interface  toughness  on  the  basis  of  the  load- 
displacement  curve,  the  size  of  the  indent,  and  the  interface  crack  length. 
Preliminary  results,  shown  in  Table  2,  indicate  that  the  interface  toughness  for  a  4  pm 
thin  film  deposited  on  a  polished  Ti-6AI-4V  surface  are  about  2.1  to  3.7  J/m2,  which 
are  in  agreement  with  those  obtained  using  the  brale  indentation  technique.  Cross- 
sectional  indentation  tests  for  the  1  pm  film  and  other  surfaces  are  in  progress. 

Table  2.  A  Summary  of  Young’s  Modulus,  Hardness,  and  Yield  Strength,  and  of 
FeCo  Thin  Films  Determined  by  Nanoindentation  Together  with 
Interface  Toughness  (Gj)  Determined  by  Cross-sectional  Nanoindentation 


Test  Sample 

Young’s  Modulus  of 
Film  (GPa) 

Hardness  of  Film 
(GPa) 

Yield  Strength  of  Film 
(MPa) 

G„  J/m2 

4.1  pm  FeCo  Film 
on  polished  Ti-6-4 

167.7 

10.7 

3570 

2.08  ±  1.57 

3.73  ±2.1 

4.1  pm  FeCo  Film 
on  unpolished  Ti-6-4 

164.8 

10.4 

3470 

In  Progress 

0.9  pm  FeCo  Film 
on  polished  Ti-6-4 

201.4 

8.5 

2830 

In  Progress 

0.9  pm  FeCo  Film 
on  unpolished  Ti-6-4 

206.1 

8.7 

2900 

In  Progress 

3.5.2  Debonding  Stress  for  FeCo/Ti-6AI-4V 

The  critical  stress  at  the  onset  of  interface  decohesion  in  FeCo/Ti-6AI-4V  was 
determined  as  s  a  function  of  thin  film  thickness  by  monotonic  loading  via  three-point 
bending.  The  interface  debonding  stresses  by  three-point  bending  are  summarized  in 
Table  1  together  with  the  interface  toughness  determined  by  indentation.  As  shown  in 
Table  1,  FeCo  thin  films  with  thickness  greater  than  4  pm  debonded  at  stresses  below 
600  MPa.  The  debonding  stress  decreases  with  increasing  film  thickness.  In 
comparison,  FeCo  thin  films  with  thickness  in  the  1  -  3  pm  range  remained  intact  and 
did  not  debond  at  a  tensile  stress  as  high  as  758  MPa,  implying  that  the  debonding 
stress  was  greater  than  758  MPa  for  the  FeCo  films  with  a  thickness  of  1  -  3  pm. 

In  addition  to  film  thickness,  the  debonding  stress  of  FeCo  thin  film  was  also 
determined  as  a  function  of  fatigue  cycles.  Figure  3.7(a)  and  (b)  summarize  the 
conditions  of  the  thin  film  for  various  peak  stresses  and  fatigue  cycles  for  4  pm  thin 
films  with  tensile  and  compressive  residual  stresses,  respectively.  In  Figure  3.7(a), 
the  FeCo  thin  film  remained  intact  at  stresses  below  440  MPa,  but  debonded  at 
stresses  above  440-460  MPa.  In  Figure  3.7(b),  the  thin  film  was  intact  at  stresses 
below  380  MPa,  but  started  to  form  cracks  at  stresses  above  380  MPa  and  become 
debonded  at  stresses  above  500  MPa.  For  both  cases,  the  debonding  stresses  were 
essentially  independent  of  the  number  of  fatigue  cycles.  On  the  other  hand,  the  thin 
film  with  a  compressive  residual  stress  exhibited  a  somewhat  higher  debonding  stress 
(440  MPa,  Figure  3.7(a))  than  the  one  with  a  tensile  residual  stress  (380  MPa,  Figure 
3.7(b)). 
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Figure  3.7.  Debonding  stress  values  of  4pm  thick  FeCo  thin  films  bonded  to  Ti- 
6AI-4V  substrate  as  a  function  of  the  number  of  fatigue  cycles:  (a) 
with  compressive  residual  stresses  and  no  surface  pre-cleaning,  and 
(b)  with  tensile  residual  stresses  and  surface  pre-cleaning. 


A  summary  of  the  debonding  stresses  for  FeCo  thin  films  on  Ti-6al-4V  is  presented  as 
a  function  of  thin  thickness  in  Figure  3.8,  which  indicates  that  the  debonding  stress 
increases  with  decreasing  FeCo  thickness  for  thickness  less  than  4  pm.  For  film 
thickness  less  than  3  pm,  the  thin  films  did  not  debond  at  748  MPa.  For  thin  films 
deposited  on  surfaces  without  ion-cleaning,  the  interface  toughness  was  about  2  J/m2, 
while  Gj  was  about  4  J/m2  for  thin  films  deposited  on  ion-cleaned  surfaces.  The 
corresponding  debonding  stresses  computed  based  on  a  critical  energy  release  rate  are 
in  reasonable  agreement  with  the  experimental  data,  as  shown  in  Figure  3.8.  The 
lower  interface  toughness  (2  J/m2)  gives  debonding  stresses  that  are  comparable  to  the 
lower  bound  data  of  the  thicker  thin  films  (  >  4.2  pm),  while  the  higher  interface 
toughness  (4  J/m2)  gives  debonding  stresses  that  are  comparable  to  the  upper  bound 
data  of  the  thin  films  over  the  entire  range  of  film  thickness  investigated.  The  results 
in  Figure  3.8.  suggest  that  interface  debonding  in  FeCo/Ti-6AI-4V  is  governed  by  a 
critical  elastic  energy  release  rate.  The  debonding  stress  is  seen  to  increase  with 
decreasing  film  thickness.  A  large  variability  exists  in  the  experimental  values  of  the 
debonding  stress  for  FeCo/Ti-6Al-4V  because  the  interface  toughness  varies  with 
location  and  shows  a  substantial  scatter  (25  -  40%  COV).  In  general,  there  is  a 
tendency  for  decohesion  to  occurs  at  the  location  where  the  interface  toughness  is 
lowest  assuming  the  local  stresses  in  the  neighboring  surround  are  approximately 
equal.  As  a  result,  the  observed  debonding  stress  values  often  fall  below  those 
predicted  on  the  basis  of  the  average  values  of  the  interface  toughness,  as  shown  in 
Figure  3.8. 
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Figure  3.8.  A  comparison  of  the  measured  and  predicted  debonding  stresses  for 
FeCo  thin  films  on  Ti-6AI-4V  shows  that  the  debonding  stress  in 
FeCo  thin  film  increases  with  decreasing  film  thickness  according  to 
a  critical  elastic  energy  release  rate  criterion.  Ion-cleaning  increased 
the  interface  toughness  from  2  J/m2  to  4  J/m2  and  increased  the 
debonding  stress. 

3.5.3  Fracture  Mechanisms 

For  both  monotonic  and  cyclic  loading,  the  debonding  process  appeared  to  be  driven 
by  the  peak  stress.  The  sequence  of  events  involved  in  the  film  decohesion  process 
that  have  been  identified  from  several  interrupted  monotonic  and  cyclic  tests.  The 
first  step  of  the  decohesion  process  in  FeCo/Ti-6Al-4V  was  the  formation  of  a  crack 
in  the  FeCo  thin  film.  Figure  3.9(a).  Once  a  crack  was  present  in  the  thin  film,  the 
thin  film  adjoining  the  crack  became  debonded.  The  debonded  regions  reflected  light 
differently  and  appeared  dark  in  the  optical  micrograph  shown  in  Figure  3.9  (b). 
Continued  cyclic  loading  caused  the  both  the  film  crack  and  the  debonded  interface  to 
grow  with  increasing  stresses  or  the  number  of  fatigue  cycles,  eventually  leading  to 
detachment  of  patches  of  the  thin  film  from  the  substrate,  as  shown  in  Figure  3.9  (c). 
In  all  cases,  cracks  in  FeCo  resulted  in  film  detachment  without  propagating  into  the 
Ti-6A1-4V  substrate. 
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(a) 


(b) 


(c) 

Figure  3.9.  Fracture  mechanisms  observed  in  FeCo/Ti-6Al-4V  subjected  to 
three-point-bend  fatigue  at  a  stress  ratio,  R,  of  0.1:  (a)  microcracks  in 
the  FeCo  thin  film,  (b)  interface  debonds  (dark  areas)  in  FeCo/Ti-6-4 
and  the  cooperative  growth  of  microcracks  and  interface-debonds 
with  fatigue  cycling,  and  (c)  debonded  patches  on  FeCo  thin  film. 


The  fracture  mechanics  of  crack  deflection  versus  crack  penetration  has  been 
analyzed  by  He  et  al.  [3]  who  presented  the  result  in  terms  of  a  plot  of  Gj/Gp  as  a 
function  of  the  Dundurs  parameter  [4],  a  =  (Ej  -  E2)/(Ei  +  E2)  [3],  where  Gj  is  the 
interface  toughness  or  the  critical  energy  release  of  the  interface  crack  and  Gp  is  the 
critical  energy  release  rate  for  crack  penetration  into  the  substrate;  Ei  and  E2  are  the 
Young’s  moduli  of  the  thin  film  (Ei  =  1.7  H  105  MPa)  and  substrate  (E2  =  1.7  H  105 
MPa),  respectively.  The  theoretical  crack  morphology  transition  boundary  is  shown 
in  Figure  3.10  as  the  solid  curve.  Crack  penetration  into  the  substrate  occurs  for  Gj/Gp 
ratios  above  the  solid  curve,  while  interface  cracks  remain  to  propagate  along  an 
interface  for  Gj/Gp  values  that  are  less  than  the  solid  boundary.  One  possible  measure 
of  the  Gp  is  the  fracture  toughness  of  the  substrate.  For  Ti-6A1-4V,  Kjc  is  at  least  67 
MPa/m  [5].  This  leads  to  a  Gp  in  excess  of  2.4  x  104  J/m2.  The  value  of  G,  range 
from  2-4  J/m2  for  FeCo/Ti-6Al-4V.  The  corresponding  Gj/Gp  ratio  is  8E-5  to  1.6E- 
4,  which  are  shown  as  filled  circles  in  Figure  3.10.  Another  possible  measure  of  Gp  is 
the  fatigue  crack  growth  threshold,  AKth,  of  small  cracks  or  large  cracks  at  a  high  R 
ratio.  For  Ti-6-4,  AKth  =  2  MPa/m  [5]  and  it  leads  to  Gp  =  21  J/m2  and  a  Gj/Gp  ratio 
of  0.095  -  0.19,  which  are  shown  as  filled  squares  in  Figure  3.10.  For  both  measures, 
the  ratios  of  Gj/Gp  ratios  are  below  the  boundary  for  crack  penetration  into  the 
substrate.  Thus,  the  interface  cracks  are  expected  to  propagate  along  the  interface  and 
would  not  penetrate  into  the  substrate.  On  the  other  hand,  the  Gj/Gp  ratio  for  a  fatigue 
crack  is  0.19,  which  is  not  too  far  the  critical  value  of  0.25  required  for  crack 
penetration  to  occur.  An  interface  toughness  in  excess  of  5.25  J/m2  would  move  the 
Gi/Gp  ratio  on  or  above  the  transition  boundary.  Thus,  the  interface  toughness  of 
FeCo  thin  films  are  currently  optimized  for  interface  fracture  without  fatigue  crack 
penetration  into  the  Ti-alloy  substrate.  The  finding  is  supported  by  experimental 
observations,  shown  in  Figure  3.9(c),  which  indicates  that  the  interface  cracks 
propagated  across  the  entire  width  and  length  of  the  FeCo  thin  film  without 
penetrating  into  the  substrate.  In  many  instances,  small  patches  of  the  debonded  thin 
films  peeled  away  from  the  substrate  while  still  attached  to  the  surface.  In  others 
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cases,  the  patches  became  detached  from  the  specimens  leaving  an  opening  on  the 
thin  film. 
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Figure  3.10.  Plot  of  Gj/Gp  ratio  as  function  of  the  Dundurs  parameter  a  [4J 
shows  the  interface  toughness,  Gj,  of  FeCo  thin  films  is  small 
compared  to  the  critical  energy  rate  for  crack  penetration,  Gp,  into 
Ti-6AI-4V  such  that  the  Gj/Gp  ratios  reside  in  the  region 
characterized  by  interface  debonding  exclusively  and  is  far  from 
the  theoretical  boundary  [3J  for  transitioning  to  crack  penetrating 
into  the  substrate. 

3.5.4  Fatigue  Life  Response 

The  effects  of  FeCo  thin  films  on  the  low-cycle  fatigue  life  (104  -105  cycles)  response 
of  Ti-6A1-4V  were  investigated  by  performing  three-point  bend  fatigue  of  rectangular 
plates  with  and  without  thin  films  at  a  peak  stress  of  748  MPa  and  R  =  0.1  until 
fracture.  Table  1  summarized  the  fatigue  lives  for  specimens  with  and  without  thin 
films.  For  uncoated  specimens,  the  fatigue  life  was  62.9  kcycles,  compared  to  74.6, 
66,  and  57.4  kcycle  for  1,  4,  and  7.5  pm  film  thickness,  respectively.  The  results 
indicate  that  the  fatigue  life  of  Ti-6A1-4V  with  thin  films  is  essentially  identical  to 
those  without  thin  films.  For  the  coated  specimens,  the  thin  films  appeared  to  have  a 
small  effect  on  the  fatigue  life,  with  increasing  fatigue  life  with  decreasing  film 
thickness.  Examination  of  the  fatigue  surfaces  revealed  that  fatigue  cracks  initiated  at 
the  edge  of  the  bend  specimens  for  both  coated  and  uncoated  specimens  with  similar 
fracture  surfaces.  In  all  four  cases,  fatigue  cracks  initiated  at  the  comer  edge  of  the 
specimens,  but  the  grain  where  the  fatigue  crack  first  appeared  was  not  identified 
because  of  the  absence  of  any  identifiable  markings  or  discemable  features.  For  the 
coated  specimens,  fatigue  cracks  initiated  on  the  side  edges  of  the  specimens  and  was 
not  located  at  or  beneath  the  FeCo  thin  film.  Thus,  fatigue  crack  initiation  in  the 
FeCo-coated  specimens  could  not  be  attributed  to  the  FeCo  thin  films.  The  similarity 
of  the  fatigue  surfaces  suggests  that  fatigue  cracks  initiated  in  the  Ti-6A1-4V  substrate 
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for  both  coated  and  uncoated  materials,  which  explains  the  similarity  in  the  observed 
fatigue  lives. 

The  fatigue  life  of  the  Ti-6A1-4V  plate  was  previously  characterized  by  axial  fatigue 
of  rounded  bars  as  part  of  the  AFRL  High  Cycle  Fatigue  (HCF)  Program  [6].  The 
fatigue  life  results  from  this  study  are  compared  against  those  for  the  HCF  program  in 
Figure  3.1 1,  which  shows  the  three-point  bend  fatigue  data  from  this  study  exhibited 
slightly  longer  life  than  the  axial  fatigue  life  data  from  the  HCF  program.  The 
difference  in  the  bending  and  axial  fatigue  lives  is  attributed  to  the  stress  gradient 
associated  with  three-point  bend  fatigue  and  not  from  the  presence  of  thin  films  on  the 
bend  specimens. 


CL 


1000 

900 


<2  800 

£ 

CO  700 

E 

|  600 

S 

5  500 


400 


R  =  0.1 


aa. 


104 


Experimental  Data 

FeCo/Ti-6AI-4V 

Ti-6AI-4V 

♦  7  |im  (Debonded) 

A 

HCF  Program 

□  4  fim  (Debonded) 

O  3  fim  (Bonded) 

□ 

Present  Study 

®  1  nm  (Bonded) 

A  A 


A 

*AA 


A  A 

/A  A 


A 

A 

A 


A4 

-£ 


105 


10® 


107 


Nf,  Cycles 


Figure  3.11.  Fatigue  life  of  three-point-bend  fatigue  of  Ti-6A1-4V  with  and 
without  FeCo  thin  film  compared  to  those  of  axial  fatigue 
generated  in  the  Air  Force  HCF  Program  [6J. 


3.5.5  Conclusions 

The  conclusions  reached  in  this  investigation  are  as  follows: 

1.  The  debonding  stress  of  FeCo  thin  films  on  Ti-6A1-4V  increases  with 
decreasing  film  thickness  according  to  a  critical  energy  release  rate. 

2.  The  interface  toughness  of  FeCo/Ti-6Al-4V  is  about  2-4  J/m2  for  film 
thickness  ranging  from  1  pm  to  8  pm. 

3.  The  interface  toughness  of  FeCo  thin  film  on  Ti-6A1-4V  varies  with  location 
and  interface  decohesion  tends  to  occur  at  the  location  with  the  lowest 
interface  toughness. 

4.  Ion-cleaning  of  Ti  surfaces  can  enhance  the  interface  toughness  and  adhesion 
of  FeCo  thin  film  on  Ti-6A1-4V. 
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5. 


The  sequence  of  damage  mechanisms  in  the  FeCo  thin  films  is  the  formation 
of  microcracks  in  the  thin  films,  followed  by  interface  crack  formation  and 
propagation,  and  then  detachment  of  thin  film  patches  from  the  substrate. 

6.  The  presence  of  thin  FeCo  films  on  Ti-6AI-4V  has  no  effect  on  the  fatigue  life 
of  the  Ti-6A1-4V  substrate  in  the  low-cycle  fatigue  regime. 

7.  FeCo  thin  films  are  durable  as  magnetostrictive  sensors.  The  overall 
durability  of  the  FeCo  thin  films  is  in  the  range  needed  for  practical 
application  of  the  film  as  imbedded  sensors. 

Though  not  presented,  experimental  evidence  in  this  study  has  also  led  to  the 

conclusion  that  the  presence  of  microcracks  and  interface  debonds  in  the  FeCo  thin 

films  did  not  affect  the  functionally  of  the  sensor  to  detect  strain  via  the  inverse 

magneto-elastic  effect. 
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